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ABSTRACT 


GAS-SOLID  TRANSPORT  IN  A  0.0508  m  PIPE  AT  VARIOUS 
INCLINATIONS  WITH  AND  WITHOUT  ELECTROSTATICS 


Craig  A.  Myler,  M.S. 

University  of  Pittsburgh 


The  transport  of  solid  particles  by  air  through  a  0.0508  m  pipe  was  studied  in  vertical, 
horizontal,  and  45°  orientations.  Through  control  of  the  air  humidity,  the  effects  of 
electrostatic  charging  was  observed.  Pressure  drop  and  particle  velocities  were  measured. 
Particles  used  included  79pm,  125pm,  and  450pm  glass  beads  and  128pm  Plexiglas 
beads. 

Analysis  of  particle  velocity,  pressure  drop,  pressure  drop  fluctuation,  electrostatic 
pressure  drop,  choking,  and  saltation  was  performed.  Visual  observations  of  the  flow 
patterns  and  behavior  were  made.  A  linear  stability  analysis  for  the  three  orientations 
was  performed. 
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NOMENCLATURE 


Drag  coefficient  of  single  particle  (dimensionless) 

Particle  diameter  (m) 


Tube  diameter  (m) 


Gas  friction  factor  (dimensionless) 


Particle  friction  factor  at  choking  (dimensionless) 

Solid  friction  factor  (dimensionless) 

Gravitational  acceleration  (m/s') 


Mass  of  gas  in  differential  length  of  tube  (kg) 

Mass  of  solids  in  differential  length  of  tube  (kg) 
Electrostatic  pressure  drop  (Pa) 

Total  pressure  drop  (Pa) 


Charge /particle  (C) 


XX 


Solids  loading  (kg  solids/ kg  air) 

Particle  Reynolds  number  (dimensionless) 
Terminal  Reynolds  number  (dimensionless) 

Time  (s) 

Actual  fluid  velocity  (m/s) 

Steady  state  fluid  velocity  (m/s) 

Fluid  velocity  fluctuation  (m/s) 

Actual  fluid  velocity  at  choking  (m/s) 
Superficial  gas  velocity  (m/s) 

Superficial  gas  velocity  at  saltation  (m/s) 
Superficial  gas  velocity  at  choking  (m/s) 
Particle  velocity  (m/s) 

Steady  state  particle  velocity  (m/s) 


Particle  velocity  fluctuation  (m/s) 

Slip  velocity  (m/s) 


XXI 


U  Terminal  velocity  (m/s) 

( 

V  Volume  of  a  particle  (m3) 

p 

Z  Contact  distance  (m) 

0 

W  Solid  flow  rate  (kg/s) 

S 

Greek  Letters 


e  Gas  voidage  (dimensionless) 

e  Gas  voidage  at  choking  (dimensionless) 

c 

e  Permittivity  of  free  space  (farad/m) 

0 

Pf  Fluid  density  (kg/m3) 

p  Particle  density  (kg/m3) 

p 

U  Fluid  viscosity  (kg/ms) 


_ _  Pressure  gradient  in  axial  direction  (Pa/m) 

:x 

ir  Constant  equal  to  3.14159265 

9  Angle  from  the  horizontal  (degrees) 


r  vVv'/VvV^i  v 1 


XXII 


Particle-Particle  collision  factor  (dimensionless) 


Friction  factor  (dimensionless) 


Wall  collision  factor  (dimensionless) 


1.0  INTRODUCTION 


In  the  transport  of  solid  partir'js  by  a  gas  stream,  numerous  forces  act  on  the 
system.  The  direction  and  magnitude  of  these  forces  significantly  affects  the  behaviour 
of  the  flow.  Two  such  forces  which  arise  in  gas-solid  transport,  and  not  significantly 
in  single  phase  gas  flow,  are  those  of  gravity  and  electrostatics.  The  major  problem 
with  these  two  forces  is  that  they  affect  the  solid  particles  almost  exclusively. 

Gravity  imposes  a  constant  downward  acceleration  on  the  particles.  To  overcome 
the  force  caused  by  this  acceleration,  the  gas  must  equal  and  overcome  this  force  by 
exerting  a  drag  force  on  the  particles.  When  the  drag  force  is  no  longer  sufficient  to 
overcome  the  force  of  gravity,  and  other  forces  present,  the  particles  fall,  and  transport 
ceases.  This  falling  of  the  particles  is  dependent  on  the  orientation  of  the  pipe  through 
which  they  are  being  transported.  In  vertical  upflow,  the  particles  can  fall  unimpeded 
by  the  pipe  itself,  whereas  in  horizontal  flow,  the  maximum  descent  is  the  diameter  of 
the  pipe.  At  other  orientations,  the  particles  can  continue  to  fall  while  in  contact  with 
the  pipe  wall.  The  effect  of  gravity  is  different  in  different  orientations  and  therefore 
must  be  accounted  for  in  conjunction  with  that  orientation. 

Electrostatic  forces  are  not  as  easily  defined.  The  contacting  of  the  solid  particles 
with  the  pipe  wall  causes  static  electrification.  As  the  particles  themselves  are  moving 
and  charged,  the  magnitude  and  direction  of  the  electric  field  is  difficult  to  define, 
even  for  a  homogeous  dispersion  of  the  particles.  Combined  with  the  effect  of  gravity, 
the  electrostatic  forces  can  affect  the  flow  differently  with  pipe  orientaiton. 

The  purpose  of  this  study  is  to  examine  the  affects  of  pipe  orientation  on  gas- 
solid  transport  with  and  without  electrostatics.  Additionally,  the  use  of  electrostatic 
ring  probes  for  the  measurement  of  solid  particle  velocity  in  a  0.0508  m  pipe  will  be 
assessed. 


2.0  LITERATURE  REVIEW 


Investigations  into  pneumatic  transport  have  generally  taken  two  approaches.  The 
first  is  an  attempt  to  empirically  correlate  data  from  these  systems  into  workable  ex¬ 
pressions.  The  second  approach  is  through  theoretical  modeling.  The  first  approach 
suffers  from  the  number  of  variables,  and  therefore,  the  number  of  groupings  required 
to  describe  a  given  system.  The  second  approach  is  hampered  by  a  lack  of  fundamen¬ 
tal  knowledge. 


11  BALANCE  OF  FORCES 

An  application  of  Newton's  second  law  to  the  particles  shown  in  figure  2-1 
provides  the  following  equation: 


Am  dU  /dt  = 

p  p 


ur 

D 


-  dF  sin  9 

£ 


-  dF  -  Am  /p  3P/3x  ±  F 

r  p  p  »dd 


(2-1) 


Similarly,  the  gas  must  also  have  a  balance  of  forces,  which  can  be  described  by: 


dF  =  3/4C  p  (U  -U  )Am  /  [(p  -p  )D  ] 

D  DS  t  ft  p»  P  p  1  P 


Wen  and  Yu'  suggest  a  modification  to  the  single  particle  drag  coefficient  of: 


C  (modified)  =e  C  (single  particle) 

DS  DS 


This  gives: 


dFo  >  3/4  £''  CKp(U  -IMA®  /  t(P(-p)D  ] 


The  second  term  on  the  RHS  of  equation  2-1  is  a  gravity  force.  Direct  applica¬ 
tion  of  Newton’s  law  gives: 


Parenthetical  references  placed  superior  to  the  line  of 
text  refer  to  the  bibliography. 
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which  is  the  weight  of  the  particles  in  the  tube. 

The  third  term  on  the  RHS  of  equation  2-1  is  the  force  due  to  friction.  Using  a 
friction  factor  for  the  solids,  this  force  can  be  described  as: 


dF  =  2f  U  2 Am  /D 

f  s  p  P  t 


(2-7) 


The  fourth  term  on  the  RHS  of  equation  2-1  is  the  force  due  to  the  pressure 
gradient  along  the  axis  of  the  pipe. 

Additional  forces  are  present  in  the  flow  system  and  are  often  significant  They 
include  forces  due  to  electrostatic  generation,  external  fields  of  force  (such  as  magneto 
static  fields),  lift  forces  due  to  rotation,  and  cohesive  forces.  The  electrostatic  forces 
will  be  discussed  later.  The  other  forces  are  assumed  to  be  negligible  for  this  discus¬ 
sion  and  will  be  neglected. 


Substitution  of  the  expessions  above  into  equation  2-1  and  rearrangement  yields: 


2f  U  */D  -  i/p  3P/3x 

*  P  l  n 


For  steady  flow,  and  neglecting  the  pressure  term,  this  becomes: 


3/4e^7C  p(U-U)2/[(p  -  p)D] 
os  f  r  p  p  f  p 


gsine  -  2f  U  VD  =  0 

*  r  1 


(2-9) 


As  the  voidage  can  be  expressed  as: 


e  =  1  -  4W/(P  irD  :U  ) 

»  p  t  p 


(2-10) 


equation  2-9  can  be  made  explicit  for  the  particle  velocity.  There  are.  however,  two 
terms  in  equation  2-9  which  make  the  solution  for  the  particle  velocity  difficult  They 
are  the  drag  coefficient  and  the  solid  friction  factor.  The  drag  coefficient  is  dependent 
on  the  air  flow  rate.  The  friction  factor  is  apparently  strongly  influenced  by  pipe 
orientation,  particle  diameter,  particle  and  pipe  material,  and  possibly  other  factors. 


2.2  PRESSURE  DROP 


The  pressure  drop  across  a  length  of  a  pneumatic  transport  line  is  a  very  impor¬ 
tant  design  parameter.  The  introduction  of  solids  into  a  gas  stream  causes,  for  most 
conditions,  an  increase  in  the  pressure  drop  over  that  of  air  alone.  Cases  of  pressure 
drop  reduction  to  below  that  of  air  alone  have  been  reported*-  .  however,  the  con- 
i  didons  for  this  phenomena  are  very  specialized.  For  steady  flow  conditions,  the  pres¬ 

sure  drop  can  be  obtained  by  adding  the  reduced  forms  of  equations  2-1  and  2-2  as  a 
sum  of  the  forces  in  the  system.  This  gives; 

Am  gsine  +  Amgsine  +Am2fU2/D  (2-11) 

P  t  p  s  p  l 

j 

+  Am2fU2/D  +(Am/p+Am/p)3P/3x=0 

!  t  *  f  1  P  P  *  I 

Note  that  in  formulating  equation  2-11.  the  force  due  to  drag  on  the  particles  is  equal 
and  opposite  the  force  causing  the  drag  by  the  fluid,  and  therefore,  cancels.  Also,  the 
additional  force  terms  have  been  omitted.  The  mass  terms  in  equation  2-11  can  be  ex¬ 
pressed  in  terms  of  the  voidage  as: 


Am  =  (l-e)p  V 

p  p  i 


(2-12) 


Am  =  ep  V 

g  g  < 


(2-13) 


and  equation  2-11  becomes: 
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(l-e)p  gsine+ep  gsine  (2-14) 

p  % 

+(l-e)p  2fU:/D 

p  «  p  i 

+Ep  2f  U  2/D  +  3P/ 3x=0 

g  *  f  i 


Integrating  equation  2-14  over  a  length  L,  rearranging,  and  noting  the  sign  convention 
of  the  pressure  drop,  yields: 


AP=  [(l-e)p  *ep  ]  Lgsine  (2-15) 

p  t 

>2p  f  (I-e)LU  2/D 
p  »  p  > 

+2p  f  eLU2/D 

p  g  f  i 


The  apparent  simplicity  of  equation  2-15  is  misleading  First,  the  particle  velocity  must 
be  known.  As  mentioned  before,  the  solids  friction  factor  is  apparently  influenced  by 
many  factors,  and  therefore  makes  solution  for  the  particle  velocity  from  equation  2-8 
difficult.  Measuring  the  particle  velocity  is  also  difficult  Methods  used  have  included 
the  use  of  radioactive  tracers*,  electrostatic  signal  cross-correlation5'  6'  ,  ultrasonic 
cross-correlation’,  and  Laser-Do ppler  veloctmeters.Q  All  of  these  methods  have  draw¬ 
backs  and/or  difficulties  in  measuring  the  particle  velocity.  The  additional  forces 


which  were  left  out  of  equation  2-15  are  sometimes  very  significant  The  electrostatic 
forces  which  can  be  generated  in  pneumatic  systems  are  often  very  large  and  have  a 
significant  impact  on  the  pressure  drop.  These  forces  will  be  discussed  further  on. 

2.3  PARTICLE  PATH  APPROACH 

Molerus10  has  taken  the  force  balance  on  particles  in  pneumatic  transport  a  step 
further.  He  has  considered  the  path  of  flight  of  a  particle  and  the  different  inter¬ 
actions  of  the  particle  in  different  phases  of  its  motion.  A  particle  is  considered  to 
undergo  the  following  flight  phases: 

Flight  under  the  influence  of  gravity 
Particle-Particle  Collision 
Particle-Wall  Collision 
Slide  along  the  pipr  wall 
Pressure  Gradient  effects 

From  these  possible  interactions,  equation  2-16  is  obtained. 


AP/p  +3p  C  (Re  )U  2L/(4D  p  )=  (2-16) 

p  f  DS  p  f  p  p 

(l-p  /p  )gL(sine+iJj  cose) 

f  p  r 

~i>  (p  u/p  )1,3lu  :/d 

p  r  p  g  p 

LU  :/D 

*  g  i 


where 


¥  =  friction  coefficient 

r 

¥  =  Particle-Particle  collision  coefficient 

p 

¥  =  Wall  collision  coefficient 

w 


Comparison  of  equation  2-16  with  equation  2-8  at  steady  state  a  obvious.  The  major 
difference  is  in  the  frictional  terms.  Molerus  provides  simplifications  to  equation  2-16 
for  different  flow  conditions;  however,  the  empirical  nature  of  the  friction  factors  re¬ 
quired  remains. 

2.4  THERMODYNAMIC  ANALOGY  APPROACH 

A  different  approach  to  describing  pneumatic  systems  was  initiated  by  Tuba." 
This  approach  treats  the  pneumatic  system  by  a  thermodynamic  phase  equilibrium  anal¬ 
ogy.  The  solids  flux,  fluid  flux,  and  voidage  are  used  in  the  format  of  the  Van  der 
Waals  equation  of  state.  In  this  format,  the  critical  properties  of  the  system,  and 
hence  the  constants  for  the  equation,  can  be  determined. 

2.5  VERTICAL  SYSTEMS 

2.5.1  Force  Balance  on  Vertical  System 


The  force  balance  equation  for  the  particles  in  a  verticle  section  of  pipe  is: 


3 


Am  dU  /dt  =  dF  -  dF  -  dF  -  Am  /p  3P/3x  ±  F 

p  p  D  t  1  P  p  add 


(2-17) 


For  dilute  systems,  the  pressure  term  is  normally  negligible  due  to  the  relatively  large 
particle  density  and  small  amount  of  particles  in  the  system.  Neglecting  additional 
forces  and  substituting  the  appropriate  expressions  for  the  forces,  this  expression  be¬ 


comes: 


dU  /dt=  3/4  £  C  p(U-U)7[(p  -  p  )D  ]  -g- 

p  DS  f  f  p  p  f  p 


2f  U  '/D  "  1/p  3P/8X 
5  r  i  p 


(2-18) 


At  steady  state,  and  with  values  for  the  drag  coefficient  and  solids  friction  factor,  this 
equation  can  be  solved  for  the  particle  velocity. 

The  pressure  drop  equation  is: 


AP=  [(l-e)o  *eo  1  Lg 

P  8 


-2d  f  (l-e)LU  */D 


-2d  f  eLL'  '  D 


(2-19) 


nr 


Solution  for  the  pressure  drop  requires  a  gas  friction  factor,  f ,  which  can  be  obtained 

g 

from  single  phase  correllations  such  as  the  Blasius  or  Koo  equations.  The  particle 

velocity  and  solids  friction  factor  are  not  so  easily  determined.  Table  2-1  lists  some  of 

the  expressions  available  and  the  systems  from  which  they  were  obtained.  The  Institute 

of  Gas  Technology  performed  testing  of  various  correlations  and  recommends  the 

modified  Konno-Saito  correlation. 

Table  2-1:  Various  Correlations  Available  for  Pneumatic 

Systems 


INVESTIGATOR 

SYSTEM 

RESULTS 

Konno  and  Saito12 

D  =.l  to  1.0  mm 

f  =  0.0285(gD)‘;:/U 

p 

p  =1440  to  2500  kg/m3 

p 

D  =26.5  and  46.8  mm 

t 

Vertical  and 

Horizontal 

S  l  P 

13 

f  = 

Vertical  s 

D  =6.78  and  13.5  mm 

l 

0.00515(l-e)/e3C(l-eU  /U 

Yang 

t  r 

Horizontal 

f$  =  0.02925(1— e)/ e3  x 

D=50.8  and  76.2  cm 

l 

[(1-eU/U  U/(gD)u:]  *M! 

i  p  f  t 

Leung  and  Wiles14 

Vertical 

Avg  of  results 

from  van  Swaaij, 

f  =  0.05/U 

Reddy  and  Pai, 
and  Konno-Saito 

s  P 

Van  Swaaij,  et.ai.‘? 

D=0.18  m 

l 

f  =  0.08/U 

s  P 

Reddy  and  Paf6 

D=0.10  m 

[ 

f  =  0.046/U 

D  =100  to  270  ym 

P 

s  P 

f 

a 

•3 

,»TS 


Capes  and  Nakamura 


Stemerding 


D  =0.0381  m 

l 

D  =256  to  3400 

p 

p  =0.911  to  7.7  gm/cm3 

p 

D  =0.0508  m 

t 

D  =20  to  150 

p 

p  =1.6  gm/cm3 


f  =  0.048  /U 


f  =  0.003 


* 

•v** 

.f.V 

«L 

r 

& 


Molerus 


Morikawa  and 


_  ..20 

Tsujj 


Horizontal 
D=.04  and  .01  m 


Vertical 

D  =40  mm  (acrylic) 

t 

D  =1.11  to  3.43  mm 

p 

p  =923  to  969  kg/m3 

p 


State  Diagram 
from  which  P 
can  be  read 


f  =  1.503OJ  /(gD  )W2)-‘ 831 

s  p  t 


Horizontal 
same  as  above 


f  =  0.805(U  /(gD)‘,V  “3 

S  p  t 


Inclined 
same  as  above 
(30,45.and 

60  degrees) 


Figure  from 
which  f 

S 

can  be  read 


Marcus,  etal. 


Horizontal 
D  =0.1  m 

t 

D  =30|im 

p 

p  =1500  kg/m3 

p 


U  =  U(1-0.0221D  °  3p  °  5) 

p  g  p  p 


Vertical 

D=.267  to  1.023  in 
D  =109  to  2024ym 

p 

p  =53.7  to  169  lb/ft3 

p 


U  =U  -  [(l*2fU*/D) 

P  g  5  P 


4/3(p  -p  d  e4  /(p  C  )] I,: 

P  f  P  1  DS 


Konno  and  Saito 


12 


1  1 


IGT‘ 

(Modified  Hinkle) 


_ 21 

iGr 

(Modified  Konno 
and  Saito) 


Vertical  and 
Horizontal 


u  =u  -u 

P  £  ' 


Vertical 

Multiple  systems 


T  T  _  TT/,  r\  0  92  0  3  ■°-2tn-°54\ 

U  =  U  (1-D  p  p  D  ) 

P  8  P  P  f  I 


Vertical  AP  =  2f  p  U  "L/ g  D  + 

s  r  g  ci 

ultiple  Systems 

0.057U  p  L/(gD)‘/: 

g  f  i 


+W L/U  ^  pL 

s  p  f 


2.6  HORIZONTAL  SYSTEMS 

The  force  balance  for  the  particles  in  a  horizontal  system  (0=0)  is: 


Am  dU  /dt  = 

p  p 


dF  -  dF  -  Am  /p  3P/3x  ±  F 

D  f  P  P  ,|W 


(2-20) 


Again,  neglecting  the  pressure  term  and  additional  forces,  this  becomes: 


dU  /dt  =  l/Az*' C  p  (U -U  TAm  /(p  -p  )D  -  2f  U ‘/D  -  Am  /p  3P/  3x  (2-21) 

p  DS  f  1  p  P  pfp  SP  1  PP 


The  effects  of  gravity  are  not  immediately  apparent  in  equation  2-21.  The  effects  are 


incorporated  in  the  solid  friction  factor.  The  expression  derived  by  Molerus10  does 
contain  a  term  in  the  horizontal  case  which  is  gravity  dependent  The  pressure  drop 
equation  for  horizontal  transport  is: 


aP=  2p  f  (l-e)LU  :/D  +  2pfeLUVD  (2-22) 

P  5  pi  f  8  f  1 


Some  correlations  applicable  to  horizontal  systems  are  given  in  table  2-1.  Again,  the 
system  from  which  the  correlations  were  obtained  is  a  key  factor  as  to  the  applicability 
of  the  correlation. 

t 

2.7  INCLINED  SYSTEMS 

The  most  general  form  of  the  force  balance  on  the  particles  in  pneumatic  trans¬ 
port  is  for  the  inclined  geometry.  Here,  the  gravity  effect  is  present  in  both  a  vertical 
and  horizontal  sense.  The  force  balance  and  the  pressure  drop  equations  are  given  by 
equations  2-9  and  2-15.  There  seems  to  be  a  general  lack  of  data  on  inclined  systems, 
and  therefore,  very  few  correlations  to  describe  them.  The  work  done  by  Morikawa 
and  Tsuji20  used  data  from  pipes  at  three  inclinations.  Their  results  give  a  correlation 
in  terms  of  a  figure  of  friction  factor  and  loading  versus  particle  Froude  number.  The 
work  of  Molerus'0  is  capable  of  describing  the  inclined  system,  however,  the  difficulties 
in  obtaining  the  separate  friction  factors  in  the  equation  would  be  severe. 


18  ELECTROSTATICS 

During  the  course  of  transporting  solid  particles  through  a  pipe,  theTe  occurs  a 
substantial  amount  of  contacting  between  the  particles  and  the  pipe  wall.  If  the  par¬ 
ticles  and  the  pipe  are  of  different  materials,  then  electrostatic  charging  occurs.  The 
extent  to  which  this  charging  proceeds  is  dependent  on  many  factors.  These  include  the 
condition  of  the  pipe  wall,  the  condition  of  the  solid  particles,  the  relative  humidity  of 
the  carrier  gas,  the  particle  size,  and  the  particle  velocity.  These  conditions  often  allow 
for  a  significant  amount  of  charging.  The  effect  of  this  charging  may  in  some  cases 
be  useful,  while  in  other  cases,  detrimental.  A  particular  use  of  the  charging  m 
pneumatic  transport  is  in  measuring  devices  for  solids  flow.'*'  5  The  detrimental  effects 
normally  appear  as  increases  in  the  pressure  drop  of  the  system.  Another  detrimental 
effect  is  that  of  discharge  of  static  electricity  which  can  lead  to  explosions. 

Klinzing'3  gives  the  force  attributable  to  electrostatics  on  a  single  particle  as: 


F  =  Eq  (2-23) 


with  the  resulting  force  for  a  system  of  n  particles  as: 


f  =  r  f 

s  L  1-1  e.i 


=  r  EqAm 

i  - 1 


(2-24) 


This  force  can  be  inserted  into  equation  2-1  to  determine  the  force  balance  for  a  svs- 


tem  involving  electrostatics.  A  problem  arises,  however,  when  trying  to  evaluate  the 
charge  per  particle  and  the  electric  field  strength.  Also,  the  direction  which  this  force 
is  exerted  is  not  easily  determined. 
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Ally  has  given  a  theoretical  analysis  of  this  type  system.  He  assumes  that  the 

particles  in  the  system  are,  on  an  average,  con  tamed  at  D  /4.  His  analysis  results  m 

t 

pressure  drop  due  to  electrostatics  of: 


AP  =45(l-E)VD/(16ir2e  D  6) 

E  i  Op 


(2-25) 


with  q  given  as: 


q  = 


fTk  (3X/30t  d5 

J  “  , .  1  1  A  » 

P» 


(2-26) 


Where  £  is  representative  of  each  of  the  variables  which  q  is  a  function  of.  The  first 
limitation  of  equation  2-25  is  the  complexity  which  equation  2-26  imparts.  Addition¬ 
ally,  the  assumption  of  a  uniform  concentration  of  particles  at  D  /4  is  obviously  un- 

l 

suited  to  flow  in  horizontal  and  inclined  pipes  where  a  significant  density  gradient  can 


2.9  STABILITY 


As  the  gas  velocity  is  decreased  in  a  pneumatic  system,  the  balance  of  forces  is 
maintained  by  a  decrease  in  the  system  voidage.  There  is  a  point  at  which  the  drag 
force  which  is  suspending  the  particles  becomes  insufficient  to  balance  the  forces  of 
gravity,  friction,  and  pressure.  At  this  point  the  system  can  no  longer  be  maintained  in 


a  steady  state.  The  range  of  instabilities  which  occur  at  and  near  this  point  are 


described  in  different  ways  depending  on  the  orientation  of  the  pneumatic  system.  In 


vertical  systems  it  is  known  as  choking.  In  horizontal  systems  it  is  called  saltation.  In 


all  cases,  it  is  a  difficult  situation  to  define  and  predict 


2.9.1  Choking 


The  phenomenon  of  choking  in  vertical  pneumatic  transport  is  best  described  by 
the  pressure  drop  observed  as  the  gas  velocity  is  decreased.  Figure  2-2  shows  the 


relationship. 


M  VELOCITY 


Figure  2- 2c  Pressure  Drop  vs.  Gas  Velocity 
Near  the  Choking  Point 


As  the  gas  velocity  is  decreased  from  point  A  to  point  B,  the  pressure  drop  decreases 
to  a  minimium  at  point  B.  Further  decrease  in  the  gas  velocity  from  B  to  C  shows  an 
increase  in  pressure  drop.  Point  B  is  called  the  choking  point  This  point,  as 
described  by  Yang*  ,  is  by  no  means  as  precise  as  shown  in  figure  2-2  due  to  the 
number  of  variables  which  determine  it  Yang  gives  the  following  correlation  to  deter¬ 
mine  the  choking  parameters  in  vertical  flow: 


2gD(e  ^  -  1)/U  -U)2  =  6.81  x  105(p  /p  )2: 

1  c  fc  t  f  p 


(2-27) 


2.9.2  Saltation 

The  stability  of  a  horizontal  pneumatic  system  is  somewhat  different  than  the  ver¬ 
tical  case.  As  the  gas  velocity  is  decreased,  the  particles  begin  to  separate  to  the  lower 
portion  of  the  pipe.  Further  decrease  in  the  gas  velocity  causes  some  particles  to  ac¬ 
tually  deposit  on  the  bottom  of  the  pipe.  This  deposition  is  known  as  saltation.  If  the 
gas  velocity  is  decreased  enough,  the  particles  will  eventually  fill  the  crosssection  of  the 
pipe  and  flow  will  stop. 

Jones  and  Leung26  compared  various  correlations  for  determining  the  saltation 
velocity.  They  recommend  the  Rizk  correlation  which  is  given  as: 


\  AV  A  A  -  .  A  «*, 

*  .  *_•  r*  v  « *  .*  a  \ 


rV-V-VA 


55? 


where:  D  is  in  mm. 

p 


This  correlation  was  recommended  for  it  simplicity  and  because  the  resultant  accuracy 
was  approximately  equal  to  that  of  other  more  complicated  expressions. 

2.9.3  Stability  Analysis 

Another  approach  to  describing  the  stability  of  a  pneumatic  transport  system  is 
through  the  use  of  the  basic  dynamic  equations  of  the  flow.  The  velocities  of  the  gas 
and  solids  can  be  expressed  as  a  steady  state  term  plus  a  fluctuation  from  that  steady 
state  by: 


U  =  U  *  (j  (2-29) 

P  P*  P 

u  =  u  «•  6 

I  gs  < 


If  the  fluctuating  terms  can  be  shown  to  decay,  then  the  system  is  said  to  be  stable. 
If  the  fluctuating  terms  grow,  then  the  system  is  said  to  be  unstable. 
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a.  Linear  Stability.  The  linear  stability  approach  utilizies  Taylor  Series  expansion  to 
linearize  the  non-linear  terms  in  equations  2-8  and  2-2.  Klinzing2  first  performed  this 
operation  on  a  vertical  system  without  electrostatics,  and  then  on  a  system  containing 
electrostatics.  The  result  of  this  analysis  was  the  following  second  order  differential 
equation: 


d2U  /dt2  -  X  dU  /dt  ♦  X  U  =  X  (2-30) 

P  1  P  2  p  3 

where: 


X  =  a  -  a  -  b  -  b 

112  12 


X  =  a  b  +  ab  -  ab 

2  2  1  2  2  1  2 


X  =  a  b  -  ab  -  ab 

3  10  0  1  0  0 


and 


a 

0 


a 

i 


a 

2 


3/4 e~*C  p(U  -U  )/[(p  -p)D]  -  gsine  -  2fU  2/D  -  1/p  3P/3X 

DS  f  fs  ps  p  f  P  s  ps  i  P 

6e*  7C  p  (U  -U  )/4p  D 

DS  f  Is  ps  p  p 

-4f  U  /D 

i  p<  i 


b  =  -3/4 e^C  (1-e)(U  -U  )2/4eD  "  gsinO  ~  2f  U  2/D  -  1/p  3P/9x 
0  ds  f«  p  i  i»  i  : 

b  =  -3eM1C  (1-e)(U  -U  )/4eD 

1  DS  fs  pt  p 

b  *  -2fU  /D  ' 

2  g  fs  i 

The  eigenvalues  of  equation  2-30  determine  the  stability  of  system  without  electrostatics. 
They  are  determined  from: 

m  =  [— y  ±(y  2-4y  U2]  /2  (2-31) 

1.2  1  1  2 

For  the  system  to  be  linearly  stable,  n^  and  m^  must  both  be  negative.  By  including 
the  force  due  to  electrostatics  in  a  vertical  system  given  by  equation  2-25,  the  analysis 
can  be  extended  to  systems  containing  electrostatics. 
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b.  Liapunov  Stability.  Another  approach  to  determining  the  stability  of  a  pneumatic 
transport  system  is  through  the  use  of  the  second  method  of  Liapunov.  This  method 
was  first  applied  to  pneumatic  transport  by  Joseph*1.  This  method  determines  regions 
of  stability  around  a  steady  state  value. 


m* 


3.0  EXPERIMENTAL  METHODS 


This  study  investigates  the  flow  of  spherical  solid  particles  in  a  0.0508  m  (2  inch 
nominal).  Excelon  (PVC)  pipe  at  three  different  inclinations.  Air  for  the  system  was 
supplied  by  a  7.5  horsepower  Roots  blower  followed  by  a  conditioning  system  capable 
of  adjusting  the  water  content  of  the  air.  The  solid  particles  were  introduced  into  the 
air  stream  by  means  of  a  live-bin  vibrating-screw  feeder  unit  The  gas-solid  stream 
was  transported  through  a  0.0508  m  (PVC)  pipe  which  included  a  3.05m  (10  foot)  test 
section  of  translucent  Excelon  pipe.  This  test  section  included  pressure  taps  at  each  end 
and  contained  two  aluminum  probes  used  for  determining  the  solid  particle  velocity. 
The  gas-  solids  stream  was  separated  in  a  cyclone  separator,  the  air  stream  being  passed 
to  the  atmosphere,  and  the  solids  stream  passed  to  a  storage  tank  mounted  on  a  plat¬ 
form  scale.  The  system  is  illustrated  in  figure  3-1. 


3.1  THE  TEST  SECTION 

The  main  test  section  consisted  of  a  schedule  40,  Excelon  pipe,  3.05  m  long  with 
an  internal  diameter  of  0.0508  m.  The  ends  of  the  pipe  were  fitted  with  two  inch 
nominal,  schedule  80,  PVC  flanges.  Pressure  taps  were  made  in  the  flanges  with  a 
3.175  mm  hole  through  to  the  lest  section. 

Probes  used  to  obtain  velocity  measurements  were  fitted  to  the  test  section  ensur¬ 
ing  continuity  between  the  probe  wall  and  the  inner  pipe  wall.  The  probes  were  con¬ 
structed  from  free  machining  aluminum  rod.  The  probe  construction  and  dimensions 
are  shown  in  figure  3-2.  The  first  pTObe  was  located  0.6096  m  from  the  downstream 
end  of  the  main  test  section.  The  second  probe  was  located  0.6605  m  upstream  of  the 
first  probe. 


Connersville,  Model  2406J,  Whispair  Blower,  mourned  on  a  baseplate.  The  blower  was 
powered  by  a  7.5  horsepower,  220  V,  3  phase  electric  motor  through  a  drive  belt  pulley 
system.  The  blower  operated  at  3600  RPM,  delivering  250  SCFM  of  air  at  6  psig.  As 
the  blower  is  a  positive  displacement,  rotary  lobe  type  blower,  a  pulsation  in  air 
delivery  was  encountered  with  a  frequency  of  240  Hz.  This  pulsation  was  partially 
damped  by  the  humid;t>  control  columns  and  the  control  valves  used  in  the  air  delivery 


3.2.2  Humidity  Control  and  Measurement 

Air  from  the  blower  was  passed  through  two  packed  columns  used  to  control  the 
humidity.  Each  column  was  constructed  from  0.2794  m  diameter  Plexiglas  tube  with 
1.905  cm  aluminum  end  plates.  The  first  column  was  packed  with  1.905  cm  ceramic, 
lntalox  saddles.  Water  was  delivered  to  the  center  of  this  column  through  tygon 
tubing.  An  outlet  at  the  bottom  of  this  column  allowed  water  to  be  passed  con¬ 
tinuously  through  the  column.  The  second  column  was  packed  with  1.905  cm  ceramic 
lntalox  saddles  and  Dryrite,  to  adsorb  moisture  from  the  incoming  air.  At  the  top  of 
each  column  was  a  three  inch  nominal.  PVC  Gate  Valve.  The  flow  through  each 
column  could  be  controlled  by  the  positioning  of  these  valves.  Thus,  the  total  air 
stream  humidity  could  be  adjusted  by  manually  positioning  these  two  valves.  Im¬ 
mediately  downstream  of  the  humidity  control  columns  was  a  thermometer  mounted  in 
the  pipe.  Following  the  thermometer  was  a  humidity  probe  connected  to  a 
Hydrodynamics  Hygrometer.  Model  15-3050.  By  the  proper  choice  of  sensing  elements, 
the  humidity  of  the  air  stream  was  able  to  be  determined.  The  range  of  measurable 
humidities  for  this  meter  was  12  to  90%.  Ehal  readings  from  the  hygrometer  and  the 
air  temperature  are  used  with  charts  to  determine  the  percent  relative  humidity. 

3.2.3  Air  Flowrate  Control  and  Measurement 

Air  from  the  humidity  control  columns  flowed  to  a  tee-fitting  where  a  bypass 
valve  allowed  the  control  of  air  discharged  to  the  atmosphere,  thus  controlling  the 
amount  of  air  through  the  test  section.  This  valve  was  a  thTee  inch  nominal.  PVC  gate 
valve  identical  to  the  valves  used  in  the  humidity  control  section.  Another  three  inch 
nominal,  PVC  gate  valve  was  in  line  with  the  feed  line  to  the  test  section.  This  valve 
allowed  the  flow  of  air  to  the  test  section  to  be  reduced  to  zero.  One  of  either  of 
these  valves  was  required  to  be  fully  open  during  operation  of  the  blower  to  avoid 
overpressurizing  the  blower.  Prior  to  the  discharge  end  of  the  air  delivery  system  was 
an  Elster,  Model  Q160  turbine  meter,  used  to  measure  the  flow  rate  of  air  to  the  test 
section.  This  meter  had  an  8— digit  roller  counter  which  displayed  the  volume  of  gas 
passed  through  the  meter.  By  reading  the  counter  at  the  beginning  and  end  of  a 
known  time  interval,  the  volumetric  flow  rate  of  air  was  obtained. 
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3.2.4  Piping  for  the  Air  Delivery  System 

All  piping  in  the  air  delivery  system  was  three  inch  nominal,  schedule  40,  PVC. 
Flanges  to  valves  and  the  tuTbine  meteT  were  three  inch  nominal,  schedule  80  PVC.  A 
three  inch  to  two  inch  reducer  coupling  was  attached  to  the  discharge  end  of  the  air 
delivery  system  to  accomadate  the  two  inch  test  section. 

3.3  SOLIDS  FEEDING  SYSTEM 

The  solids  feeding  system  consisted  of  a  Vibra-Screw,  Inc..  Live  Bin  Volumetric 
Screw  Feeder  mounted  on  a  steel  frame.  The  bin  capacity  of  this  feeder  was  three 
cubic  feet.  The  feeder  was  modified  to  make  use  of  an  auger  type  screw  which  al¬ 
lowed  for  higher  back  pressures  on  the  screw  feeder  than  was  possible  with  the 

manufacturers  screw.  The  auger  screw  was  mounted  in  a  chuck  which  used  teflon 
tipped  set  screws  to  secure  the  screw.  If  the  screw  was  impeded  due  to  excessive  pres¬ 
sure  or  a  blockage,  the  chuck  continued  turning  without  turning  the  screw  until  the 

blockage  was  eliminated.  The  rate  of  revolution  of  the  screw,  and  thus,  the  volumetric 

feed  rate  of  the  solids,  was  controlled  by  a  dial  setting  on  the  feeder  unit. 

The  discharge  from  the  screw  feeder  was  injected  into  the  test  section  through  a 
two  inch  nominal  y-fitting.  This  allowed  the  solids  to  be  injected  downward  into  the 
test  section,  thereby  reducing  blockages  in  the  screw  feeder. 

3.4  SOLIDS  RECOVERY  AND  WEIGHING  SYSTEM 

Solids  recovery  from  the  test  section  was  accomplished  with  a  Federal  Classifica¬ 
tions  Systems  Cyclone  Separator.  The  air  stream  from  the  cyclone  was  passed  through 
an  MSA  absolute  filter  to  the  atmosphere.  The  solids  from  the  bottom  of  the  cyclone 
fell  into  a  15  gallon,  conical  bottom  storage  tank  which  was  mounted  on  a  Circuts  and 
Systems,  Inc.,  Model  sx-501  platform  scale.  A  valve  was  located  at  the  base  of  the 
storage  tank  which  allowed  the  solid  mass  rate  to  be  determined.  With  the  valve  open. 


the  solids  from  the  storage  tank  passed  by  gravity,  back  to  the  bin  of  the  screw  feeder, 
thus  completing  the  loop. 


3.5  PRESSURE  MEASUREMENT 

The  pressure  drop  through  the  test  section  was  measured  by  the  use  of  a  Viatran 
Model  215  pressure  transducer.  The  range  of  this  transducer  was  0  to  5  inches  of 
water.  Pressure  ups  were  located  at  a  ten  foot  interval  of  the  test  section.  The  out¬ 
put  signal  from  the  transducer  was  connected  to  a  Hewlett-Packard  Model  7702B  Strip- 
chart  Recorder  w..ich  provided  a  graphical  output  of  the  pressure  response. 


3.6  PARTICLE  VELOCITY  MEASUREMENT 

Measurement  of  particle  velocity  in  the  test  section  was  accomplished  using 
electrosutic  ring  probes.5  This  method  involves  the  use  of  cross-correlation  techniques 
to  determine  time  of  flight  through  a  known  distance.  The  probes  used  were  con¬ 
structed  from  three  inch  diameter,  free-machining,  aluminum  rod.  The  conuct  length 
of  each  probe  was  two  inches.  The  probe  separation  was  2.167  feet. 

Signals  from  the  electrosutic  probes  were  processed  through  two  Keithly  610C 
Electrometers.  The  output  from  the  electrometers  was  recorded  on  magnetic  upe.  This 
provided  a  simultaneous  record  of  each  pTobe  which  could  be  processed  by  computer  to 
determine  the  particle  velocity. 

The  computer  used  for  processing  the  dau  was  a  Digiul  Equipment  Corporation 
MNC/Declab-23  system.  This  system  conuined  an  analog  to  digiul  converter  through 
which  recorded  signals  were  converted  to  digiul  signal  files.  The  digiul  signal  files 
were  then  correlated  using  a  FORTRAN  computer  program.  The  computer  programs 
used  for  dau  input  to  the  computer  and  for  determining  the  particle  velocity  are  in¬ 
cluded  in  Appendix  C. 


3.7  EXPERIMENTAL  PROCEDURE 


The  experimental  procedure  involved  three  overall  steps:  Start-up,  Experimental 
Data  Acquisition,  and  Shut-down.  The  start-up  and  shut-down  procedures  were  com¬ 
pleted  at  the  beginning  and  end  of  each  operating  period.  The  steps  taken  for  ex¬ 
perimental  data  aquisition  were  completed  for  each  condition.  Refer  to  figure  3-1  for 
equipment  numbers. 

3.7.1  Start-up  Procedures 

The  following  steps  were  completed  at  the  beginning  of  each  experimental  day: 

-  Check  equipment  for  breakage,  stray  objects,  etc. 


-  Electronic  equipment  warm-up  (approximately  15  minutes) 


-  Check  pressure  transducer  calibration 


-  Install  desired  humidity  sensor 


-  Check  valves 

•  Valves  4,5,  and  6  should  be  fully  open  and  valve  16  should  be  fully- 
closed 


-  Blower  start-up 

•  Caution  should  be  made  when  starting  the  blower  that  proper  ear 
protection  is  observed.  The  noise  level  of  the  bloweT  exceeds  80 


decibels  at  5  feet. 


-  Humidity  adjustment 

*  The  humidity  meter  is  turned  to  on.  For  lower  humidity  runs,  valve  5 
is  slowly  closed  until  the  desired  humidity  is  achieved.  For  high 
humidity,  the  water  inlet  line  is  connected  to  the  water  supply  and  the 
needle  valve  is  adjusted  for  proper  water  flow.  Valve  4  is  then  slowly 
closed  until  the  desired  humidity  is  achieved. 

*  Care  must  be  taken  not  to  allow  water  to  become  entrained  in  the  air 
stream  as  damage  to  the  humidity  probe  could  result 

3.7.2  Experimental  Data  Acquisition 

Data  acquisition  was  accomplished  using  pre-printed  data  acquisition  forms.  A 
sample  form  is  shown  in  appendix  .  Preset  conditions  for  each  run  were  first  made. 
These  included  the  test  section  orientation,  humidity,  and  particle  size.  The  air  flow 
rate  was  approximated  by  the  use  of  a  manometer,  which  read  the  air  pressure  at  the 
outlet  of  the  air  delivery  system.  The  following  steps  were  taken  for  each  run: 

1.  The  date,  time,  tape  number,  and  panicle  information  were  recorded. 

2.  Valve  16  was  slowly  opened  completely. 

3.  Valve  6  was  closed  until  the  desired  air  flow  rate  was  achieved. 

4.  The  temperature  and  humidity  were  noted. 

5.  The  solids  feeder  was  turned  to  on  and  the  dial  reading  recorded. 
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6.  The  turbine  meter  counter  reading  was  recorded  and  time  was  begun  on  a 
stopwatch 

7.  The  scale  was  reset  to  zero  if  required  and  the  valve  located  at  the  bottom 
of  the  storage  tank  was  closed.  The  time  on  the  stopwatch  was  noted  for 
the  time  of  solids  weighing. 

8.  The  tape  counter  number  on  the  tape  recorder  was  noted  and  the  recorder 
started.  After  a  sufficient  time  period  (at  least  8  seconds),  the  tape  recorder 
was  stopped  and  the  counter  number  again  noted.  The  tape  recorder  was 
then  advanced  to  provide  spacing  between  signals  for  each  run. 

9.  Visual  observations  were  then  made  through  the  glass  viewing  section. 

10.  The  scale  reading  was  taken  and  the  time  period  from  the  stopwatch  was 
recorded.  The  valve  located  at  the  base  of  the  storage  tank  was  opened. 

11.  The  turbine  meter  counter  reading  was  taken  and  the  time  period  from  the 
stopwatch  recorded. 

12.  Humidity  and  Temperature  were  again  noted. 

13.  The  ending  time  of  the  experimental  run  was  then  recorded. 


14.  The  solids  feeder  was  turned  to  off  and  valve  6  fully  opened. 


3.7.3  Shutdown  Procedure 

There  were  no  unusual  shutdown  procedures  other  than  shutting  off  of  the  equi¬ 
pment,  except  in  the  case  of  high  humidity  experiments.  Fot  these  cases,  the 
humidification  column  had  to  be  purged  to  prevent  water  from  leaking  back  into  the 
blower.  This  was  accomplished  by  stopping  the  water  supply  to  column  3  and  allowing 
the  blower  to  run  with  valve  4  fully  closed,  and  valve  5  fully  open.  When  the  water 
in  the  column  was  eliminated,  valves  4  and  5  were  opened  fully  and  the  blower 
stopped. 


4.0  RESULTS  AND  ANALYSIS 


4.1  PARTICLE  SIZE 

Experiments  were  conducted  using  four  particle  sizes.  Three  of  these  were  glass 
beads  while  the  fourth  was  Plexiglas.  A  particle  size  analysis  was  conducted  on  these 
particles  prior  to  the  experiments  and  after  completion  of  the  experiments.  Table  4-1 
shows  the  results  of  the  analysis  and  the  particle  size  used.  The  results  of  the  analysis 
are  included  in  Appendix  . 


4.2  FLOW  BEHAVIOR  AND  PATTERNS 

The  three  different  test  section  orientations  showed  a  marked  difference  as  to  the 
flow  behavior  and  flow  patterns  of  the  different  particles.  There  were  some  general 
aspects  common  to  all  particles  in  a  given  orientation.  In  the  vertical  pipe,  as  the  gas 
velocity  was  lowered,  the  particles  could  be  seen  to  deviate  from  streamlined  flow  lines. 
Pulsations  occurred  in  the  system  where  denser  slugs  of  material  moved  through  the 
pipe  with  lengths  of  lesser  particle  density  between  them. 

In  the  horizontal  test  section,  decrease  in  the  gas  velocity  caused  a  radial  separa¬ 
tion  in  the  pipe.  The  bottom  of  the  pipe  had  a  higher  particle  density  than  the  top. 
even  before  saltation  occurred. 

As  the  gas  velocity  was  decreased  in  the  inclined  orientation,  a  reverse  flow  be¬ 
havior  was  observed.  Particles  formed  retrograde  dunes  which  slid  downwards  and  were 
either  reduced  by  entrainment,  or  were  destroyed  by  partially  plugging  the  entrance  to 
the  test  section  where  they  were  redispersed  into  the  gas  stream. 


tv1 


1 


Table  4-1:  Particle  Size  Analysis 


Weight  Mean  Particle  Diameter  (ym) 


791  (Glass)  90.8  67.8  (avg=79.3) 


125.0  (Glass)  133.8  100.7  (avg=117.3) 


446.3  (Glass)  447.8  343.1  (avg=395.5) 


126.7  100.9  (avg=113.8) 


(Plexiglas) 


a.  Vertical  Systems.  The  79ym  glass  beads  tended  to  pulse  through  the  system  even  at 
higher  velocities.  These  pulsations  could  be  termed  slugging;  however,  the  density  of 
the  slugs  was  only  slightly  greater  than  that  of  the  flow  between  them  and  they  oc- 
cured  at  higher  frequencies.  Definite  slugging  was  observed  at  lower  gas  velocities.  An 
additional  observation  at  lower  velocities  was  an  internal  radial  motion  of  the  particles. 
This  motion  became  more  pronounced  as  the  gas  velocity  was  further  decreased.  This 
radial  motion  set  in  at  much  higher  gas  velocities  for  the  systems  at  high  humidity.  It 
appeared  that  the  presence  of  electrostatics  in  the  system  stabilized  this  motion. 


The  125ym  glass  beads  acted  much  the  same  as  the  79ym  glass  beads  except  that 
definite  slugging  occurred  at  higher  gas  velocities.  Again,  the  radial  motions  were 
damped  by  the  presence  of  electrostatics,  but  the  formation  of  dense  slugs  occurred  at 
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higher  gas  velocities.  Wall  interactions  were  also  more  pronounced  with  the  125vim 
glass  beads.  At  low  gas  velocities,  particles  hit  the  walls  and  then  rebounded  into  the 
main  stream  of  flow  with  increasingly  lower  tragectories.  At  the  lowest  gas  velocities, 

some  particles  hit  the  wall  and  traveled  downward  in  the  pipe. 

Wall  interactions  were  most  pronounced  with  the  450ym  glass  beads.  At  lower 
velocities,  particles  in-between  slugs,  could  be  seen  to  traverse  the  width  of  the  pipe 
from  wall  collision  to  wall  collision.  Again,  the  presence  of  electrostatics  seemed  to 
dampen  the  internal  motion  of  the  flow.  The  plugging  of  the  pipe  with  these  particles 

was  a  very  rapid  phenomenon.  The  smaller  particles  seemed  to  drop  out  slowly  in 

comparison. 

b.  Horizontal  Systems.  The  79ym  glass  beads  were  the  most  susceptible  to  saltation  in 
the  horizontal  pipe.  As  gas  velocities  were  decreased,  the  flow  separated  with  the  bot¬ 
tom  of  the  pipe  having  a  higher  solids  concentration  than  the  top  (See  Figure  4-1). 
Further  decrease  in  gas  velocity  caused  some  of  the  particles  to  salt  out  on  the  bottom 
of  the  pipe  forming  blunt  nosed  islands.  A  further  decrease  in  gas  velocity  caused 
more  saltation  until  the  bottom  of  the  pipe  was  covered  with  a  layer  of  particles.  At 
the  onset  of  this  condition,  particles  began  to  salt  out  in  the  vertical  section  and  even¬ 
tually  plugged  the  vertical  pipe.  The  presence  of  electrostatics  caused  more  pulsing  in 
the  horizontal  section. 

The  I25ym  glass  beads  followed  similar  behavior  to  the  79ym  glass  beads  except 
when  saltation  caused  the  bottom  of  the  pipe  to  be  covered.  In  this  case,  the  small 
disturbances  seen  on  the  top  of  this  layer  were  not  observed.  The  layer  was  very 
smooth  and  particles  could  be  seen  to  be  lifted  from  this  layer  into  the  stream  flowing 
above. 

The  450ym  glass  beads  did  not  salt  out  in  the  horizontal  orientation.  Although 
flow  separated  to  the  bottom  of  the  pipe,  the  particles  never  actually  stayed  on  the 
pipe  wall.  There  was  much  more  activity  along  the  particle  tragectories  and  more  wall 
collisions  were  observed. 


Figure  4-1:  Flow  Patterns  in  Horizontal  Flow 


c.  Inclined  Systems.  The  79ym  glass  beads  tended  to  begin  salting  out  of  the 
horizontal  section  before  any  deposition  occured  in  the  inclined  section.  When  deposi¬ 
tion  did  begin,  it  began  very  near  the  entrance  to  the  inclined  section.  Eventually, 
retrograde  dunes  formed  and  traveled  downward  (See  Figure  4-2). 

The  125ym  glass  beads  formed  the  same  type  retrograde  dunes  as  the  79ym  glass 
beads,  except  that  the  horizontal  section  preceeding  the  45°  test  section  salted  out  prior 


Figure  4-2:  Retrograde  Dunes  in  Inclined  Flow 


to  formation  of  the  retrograde  dunes  for  high  humidity  cases,  and  afieT  formation  foT 
the  lower  humidities. 

The  450ym  glass  beads  formed  the  retrograde  dunes  very  quickly  and  their  size 
and  velocity  was  much  higher  than  the  smaller  particles.  These  dunes  began  as  a  thin 
layer  on  the  bottom  of  the  pipe  from  which  particles  were  entrained.  This  layer  grew 
as  deposition  incrersed  until  a  dune  of  sufficient  size  formed  and  slid  down  the  pipe. 
The  formation  of  the  retrograde  dunes  in  all  systems  concluded  with  the  plugging  of 
the  vertical  section.  The  smaller  particle  sizes  took  much  longer  for  this  to  occur  than 
the  450ym  glass  beads. 


4.3  PARTICLE  VELOCITY  ANALYSIS 


The  particle  velocity  was  obtained  from  cross-correlation  of  the  signals  from  two 
electrostatic  ring  probes.  This  data  has  been  plotted  as  particle  velocity  vs.  gas  velocity. 
The  best  fit  of  the  data  to  a  straight  line  was  obtained  by  linear  regression.  Com¬ 
parison  is  made  between  the  experimental  particle  velocity,  the  expression  U  -U .  and 
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the  IGT  recommended  correlation.  The  particle  velocity  obtained  for  the  vertical  sys¬ 
tems  shows  fair  agreement  with  the  correlations.  Figure  4-3  shows  the  particle  velocity 
for  125ym  glass  beads  vs.  gas  velocity.  At  lower  gas  velocities,  agreement  with  the 
correlations  is  within  10%.  As  the  gas  velocity  increases,  the  experimental  values  tend 
to  be  less  than  that  from  either  correlation  with  a  mean  deviation  of  approximately 
20%.  For  the  79ym  glass  beads.  I25ym  glass  beads  and  128  ym  Plexiglas  beads,  the 
experimental  particle  velocity  was  less  than  both  correlations  with  the  expression  U  -U 
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giving  only  slightly  better  results.  For  the  450ym  glass  beads,  the  experimental  particle 
velocities  were  between  the  two  correlations  with  the  IGT  over  predicting.  The  ad¬ 
ditional  particle  velocity  figures  are  included  in  figures  A-l,  A-2  and  A-3. 

The  horizontal  systems  showed  much  the  same  results  as  the  vertical  systems,  ex¬ 
cept  that  the  spread  of  the  data  was  much  greater.  The  analysis  for  the  125yim  glass 
beads  is  shown  in  figure  4-4.  The  mean  deviation  from  the  correlations  has  increased 
from  20%  to  47%.  The  additional  horizontal  systems  are  shown  in  figures  A-4  and 
A-5. 


The  inclined  systems  also  show  the  same  type  deviation  from  the  two  correlations. 
The  analysis  for  the  125 ym  glass  beads  is  shown  in  figure  4-5.  The  additional  inclined 
systems  are  shown  in  figures  A-6  and  A -7. 
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4.4  PRESSURE  DROP 


Various  pressure  drop  relationships  were  investigated  to  determine  the  effect  of 
pipe  orientation  and  electrostatics  on  the  flow.  These  included  the  pressure  drop  vs. 
gas  velocity  relationship,  the  pressure  drop  fluctuation  vs.  gas  velocity  relationship,  and 
the  electrostatic  pressure  drop  vs.  water  to  air  mass  ratio. 

4.4.1  Pressure  Drop  vs.  Gas  Velocity 

As  discussed  in  section  2.9,  the  pressure  drop  vs.  gas  velocity  curve  graphically 
represents  the  nature  of  a  pneumatic  system  in  terms  of  stability.  It  also  provides  the 
criterion  for  optimum  operation.  The  presence  of  electrostatics  in  the  system  affects 
the  pressure  dTop  in  different  ways.  For  this  study,  operation  at  high  humidity  was  as¬ 
sumed  to  eliminate  the  electrostatic  effects. 

a.  Vertical  Systems.  The  pressure  drop  vs.  gas  velocities  for  the  vertical  systems 
studied  are  shown  in  figures  4-6,  A-8,  A-9,  and  A-10.  Particular  attention  was  paid  to 
these  systems  during  the  experiments  to  obtain  data  near  the  choking  point.  In  most 
cases,  the  system  was  allowed  to  completely  plug  the  pipe.  The  79pm  glass  beads 
remained  in  steady  flow  at  much  lower  gas  velocities  than  the  larger  glass  beads  as  ex¬ 
pected.  The  instabilities  near  the  choking  point  were  not  as  dramatic  as  the  125pm 
glass  beads.  By  visual  observation,  a  great  deal  of  instability  in  the  system  was  seen, 
however,  the  pressure  drop  readings  did  not  indicate  the  dramatic  increase  in  pressure 
that  was  expected.  The  effects  of  electrostatics  in  this  system  were  also  not  as  ex¬ 
pected.  The  smaller  particle  size  should  have  caused  the  largest  increase  in  pressure 
drop  of  all  particles  studied.  The  results,  however,  show  only  a  very  small  increase  in 
pressure  drop  due  to  electrostatics. 

The  125pm  glass  beads  showed  a  unique  behavior  near  the  choking  point.  Figure 
4-6  shows  an  increasing  maximum  for  the  high  humidity  cases.  This  maximum  actually 
exceeds  the  pressure  drop  for  the  low  humidity  experiments  at  the  highest  mass  flow 
rate.  A  similar  phenomenon  was  found  in  a  0.0254  m  pipe  by  Zaltash”.  This 
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Figure  4-6: 


Pressure  Drop  vs.  Gas  Velocity  for  125y 
Glass  Beads  in  the  Vertical  Orientation 
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phenomenon  was  attributed  to  a  particle  size  effect.  Since  the  Plexiglas  beads  had  a 
partical  size  of  128ym,  and  this  behavior  was  not  observed  for  them,  a  density  factor 
should  be  included.  For  Zaltash’s  study,  the  group  (D  /D)(p  /p)  was  equal  to  6.06. 

p  i  p  f 

In  this  study,  the  same  group  equals  4.79,  suggesting  that  this  group  could  describe  a 
critical  particle  to  pipe  combination.  This  type  of  phenomenon  has  been  disclaimed', 
however,  most  studies  have  not  explored  the  region  of  choking  that  closely. 

The  450ym  glass  beads  showed  results  much  more  as  expected.  The  range  of  un¬ 
stable  behavior  was  not  as  broad  as  with  the  smaller  particles  and  instability  leading  to 
the  plugging  of  the  pipe  occured  very  rapidly.  For  low  mass  flow  rates,  the  electros¬ 
tatic  effect  was  negligible.  As  the  mass  flow  rate  was  increased,  the  electrostatic  con¬ 
tribution  to  the  pressure  drop  increased. 

The  128ym  plexiglas  particles  showed  very  little  electrostatic  effects.  This  is  due 
to  the  nature  of  the  particle  and  pipe  materials  having  very  similar  dielectric  constants. 

b.  Horizontal  Systems.  Three  particles  were  studied  in  the  horizontal  section.  They 
were  the  79ym,  125ym  and  450ym  glass  beads.  The  Plexiglas  particles  were  not  studied 
due  to  concern  oveT  the  explosive  hazard  they  presented.  The  range  of  gas  velocities 
was  less  than  that  of  the  vertical  systems  as  the  increased  total  line  pressure  prevented 
consistant  bead  feeder  operation  at  higher  velocities.  As  the  horizontal  lest  section  was 
preceded  by  a  2.5m  vertical  section,  the  complete  plugging  of  the  horizontal  section  w’as 
not  obtained.  The  pressure  drop  vs.  gas  velocity  data  is  presented  in  figures  4-7,  A-ll, 
and  A-12.  The  79ym  glass  beads  were  the  most  susceptible  to  saltation  due  to  the  sal¬ 
tation  gas  velocities  being  greater  than  the  choking  velocities.  The  pressure  drop  vs.  gas 
velocity  curves  are  less  smooth  than  the  larger  particles.  The  presence  of  electrostatics 
m  the  low  humidity  cases  tended  to  dampen  the  fluctuations.  The  effect  of  the 
electrostatics  on  the  magnitude  of  the  pressure  drop  was  seen  to  vary  due  to  the  fluc¬ 
tuations  present  in  the  high  humidity  cases. 

The  pressure  drop  vs.  gas  velocity  data  for  the  I25ym  glass  beads  is  shown  in 
figure  4-7.  Again,  the  presence  of  electrostatics  dampened  the  fluctuations  in  the 
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curves.  The  effect  of  the  electrostatics  on  the  magnitude  of  the  pressure  drop  is  an 
overall  decrease.  This  is  opposite  to  the  effect  found  in  the  vertical  cases.(see  Figure 
4-6) 

The  450pm  glass  beads  showed  a  much  more  stable  pressure  drop  response.  This 
is  due  to  the  450pm  glass  beads  not  salting  out  at  all  prior  to  complete  plugging  of  the 
vertical  section.  The  effect  of  electrostatics  was  found  to  increase  the  pressure  drop, 
although,  this  effect  was  not  as  prominent  as  the  decrease  found  for  the  125pm  glass 
beads. 


c.  Inclined  Systems.  The  particles,  mass  rates,  and  humidities  studied  in  the  inclined 
section  were  the  same  as  those  for  the  horizontal  section.  The  inclined  test  section  was 
placed  at  the  same  point  as  the  horizontal  section.  The  test  section  was  inclined  at 
45°.  Due  to  space  limitations,  no  entrance  length  after  the  horizontal  section  was  in¬ 
cluded.  The  resultant  effect  of  the  45°  bend  on  the  test  section  was  not  included  in 
the  analysis.  The  results  are  presented  in  figures  A-13,  4-8  and  A-14. 


The  79pm  glass  beads  showed  the  greatest  fluctuation  in  the  pressure  drop  vs.  gas 
velocity  curves.  Electrostatic  effects  caused  increased  fluctuations.  In  the  higher  mass 
flow  rate  cases,  the  lower  humidity  pressure  drop  was  greater  than  that  of  the  high 
humidity.  This  is  a  reversal  of  the  pattern  seen  in  the  horizontal  systems. 

The  pressure  drop  vs.  gas  velocity  data  for  the  125pm  glass  beads  is  shown  in 
figure  4-8.  Again,  a  reverse  of  the  horizontal  case  was  observed  with  electrostatics 
causing  an  increase  in  the  pressure  drop  and  fluctuations. 

The  40ym  glass  beads  had  the  least  fluctuations  in  the  pressure  drop  vs  gas 
velocity  curves.  Also,  the  presence  of  electrostatics  caused  an  increase  in  pressure  drop 
for  both  mass  flow  rates. 


d.  Combined  Svstem  Analysis.  Each  system  was  compared  for  the  effect  of  pipe 
orientation  on  pressure  drop  vs.  gas  velocity.  This  analysis  shows  how  the  pressure 
drop  for  different  particle  sizes  varied  with  pipe  orientation  as  well  as  electrostatic  ef¬ 
fects.  This  combined  data  is  shown  in  figures  A-15  through  A-26. 

For  the  79pm  glass  beads,  the  horizontal  orientation  dominated  the  vertical  case, 
except  for  the  condition  of  high  mass  flow  rates  and  high  humidity.  The  effect  of 
electrostatics  dampened  fluctuations  in  the  horizontal  case  with  little  effect  in  the  ver¬ 
tical  cases.  This  can  be  explained  by  a  charge  distribution  effect,  where  in  the  vertical 
case,  clustering  is  inhibited  by  the  charge  on  the  particles  but  pressure  drop  is  increased 
due  to  the  additional  electrostatic  force.  In  the  horizontal  cases,  the  charge  on  the 
particles  tends  to  inhibit  the  increase  in  solid  concentration  to  the  bottom  of  the  pipe. 
The  45°  inclined  pipe  showed  an  increase  in  pressure  drop  over  both  the  horizontal 
and  vertical  cases  for  the  low  humidity  cases.  For  high  humidities,  the  pressure  drops 
were  much  closer  to  those  for  the  horizontal  and  vertical  cases.  This  is  apparently  due 
to  a  combination  of  the  repulsive  forces  of  the  particles  now  having  an  axial  com¬ 
ponent  and  the  electrostatic  force  terms. 

The  125pm  glass  beads  followed  the  same  trend  as  the  79pm  glass  beads,  except 
that  now  the  separation  between  vertical  and  horizontal  cases  is  not  as  great.  The 
dominant  orientation  is  not  as  clearly  delineated,  which  was  borne  out  by  observation 
when  saltation  and  choking  occurred  simultaneously. 

The  450pm  glass  beads  also  follow  the  same  pattern,  except  that  the  vertical  sys¬ 
tem  now  dominates.  This  effect  was  observed  by  the  vertical  section  choking  above  the 
saltation  velocity  of  the  particles. 
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4.4.2  Comparison  to  Correlations 

The  experimental  pressure  drop  was  compared  to  two  correlations  for  the  vertical 
and  horizontal  orientations.  The  Konno  and  Saito1'  correlation  was  used  for  both 
orientations.  The  correlations  by  Yang’3  were  also  used,  one  for  the  vertical  orientation 
and  another  for  the  horizontal  orientation.  The  results  of  this  analysis  have  been 
plotted  as  APexperimental/APcalculated  vs.  gas  velocity  and  are  shown  in  Figures  A-40 
through  A-56.  Absolute  mean  errors  and  standard  deviations  were  also  calculated  and 
are  included  in  Tables  B-l  through  B-7. 

In  the  vertical  systems,  the  correlation  of  Konno  and  Saito  had  an  absolute  mean 
error  range  of  11.0%  to  88.1%.  The  correlation  of  Yang  had  an  error  range  of  11.7% 
to  44.3%.  It  is  important  to  note  that  for  most  conditions,  the  values  calculated  using 
the  correlation  of  Yang,  were  less  than  the  experimental  values,  while  the  correlation  of 
Konno  and  Saito  usually  predicted  pressure  drops  above  the  experimental  values. 

For  the  horizontal  systems,  the  correlation  of  Konno  and  Saito  gave  an  error 
range  from  10.2%  to  55.9%.  The  correlation  of  Yang  had  an  error  range  from  6.0%  to 
36.7%.  For  most  conditions,  the  correlation  of  Konno  and  Saito  was  found  to  under¬ 
predict  the  pressure  drop.  Such  a  generalization  could  not  be  made  for  the  results  ob¬ 
tained  using  the  correlation  of  Yang. 

4.4.3  Pressure  Drop  Fluctuation 

Upon  introduction  of  the  solids  to  the  gas  stream;  the  pressure  drop  began  to 
fluctuate.  This  fluctuation  was  found  to  depend  on  particle  size,  humidity,  pipe  orien¬ 
tation,  and  solids  mass  flow  rate. 


a.  Vertical  Systems.  The  pressure  drop  vs.  gas  velocity  for  the  79ym  glass  beads  in 
the  vertical  pipe  is  shown  in  figure  A-28.  There  was  little  correlation  between  gas 
velocity  and  pressure  drop  fluctuation  with  a  mean  fluctuation  of  about  20%.  The  ef¬ 
fect  of  electrostatics  and  solids  mass  flow  rate  was  also  indeterminate. 

The  125ym  glass  beads  (see  Figure  )  showed  an  increase  in  pressure  drop  fluctua¬ 
tion  with  decrease  in  mass  flow  rate.  This  could  be  due  to  a  decrease  in  particle- 
particle  collisions  allowing  greater  panicle  mobility  across  the  pipe.  The  presence  of 
electrostatics  increased  the  pressure  drop  fluctuations  for  all  cases.  At  the  lowest  mass 
flow  rate,  a  maximum  fluctuation  was  found  as  the  gas  velocity  was  decreased.  Further 
decrease  in  gas  velocity  caused  a  decrease  in  fluctuation  intensity.  This  was  also  the 
case  for  the  next  highest  mass  flow  rate,  except  the  maximum  was  not  as  well  defined. 
At  the  highest  solids  mass  flow  rate,  this  maximum  was  not  present 

The  data  for  the  450ym  glass  beads  is  shown  in  figure  A-29.  Again,  a  decrease 
in  fluctuations  was  found  for  increasing  solids  mass  flow  rate.  For  all  cases,  fluctua¬ 
tions  increased  with  decreasing  gas  velocity.  The  presence  of  electrostatics  caused  an 
increase  in  fluctuations  for  the  two  lower  mass  flow  rates  but  a  reduction  in  fluctua¬ 
tions  for  the  higher  solids  mass  flow  rate.  In  all  cases  this  difference  was  less  than 
10%. 

The  128ym  Plexiglass  beads  showed  little  fluctuation  in  pressure  drop,  (see  Figure 
A-30)  A  mean  fluctuation  of  about  5%  best  describes  this  system.  This  is  pTObably  due 
to  the  lesser  density  of  the  Plexiglas  panicles  (about  1/2  that  of  the  glass  beads)  and 
the  low  mass  flow  rates  studied. 

b.  Horizontal  Systems.  The  results  for  the  79ym  glass  beads  in  the  horizontal  orien¬ 
tation  are  shown  in  figure  A-32.  Like  the  vertical  case,  the  fluctuations  appear  to  van 
about  a  mean  value  of  20%  with  little  effect  from  gas  velocity  or  solids  mass  flow  rate 

The  results  for  the  I25ym  glass  beads  are  shown  in  Figure  A-31.  In  this  case 
there  is  an  increase  in  pressure  drop  fluctuation  with  decreasing  gas  velocity.  For  the 


higher  solids  mass  flow  rate,  the  presence  of  electrostatics  increases  the  fluctuations,  but 
only  by  about  2%.  In  the  lower  mass  flow  rate  there  is  not  appreciable  difference  be¬ 
tween  the  electrostatic  and  non-electrostatic  cases. 

The  results  for  the  450ym  glass  beads  are  shown  in  Figure  A-33.  There  is  an  in¬ 
crease  in  pressure  drop  fluctuation  with  decreasing  gas  velocity.  Solids  mass  flow  rate 
had  little  effect  on  the  pressure  drop  fluctuations  while  electrostatics  showed  a  small 
decrease  in  fluctuations. 

c.  Inclined  Systems.  The  inclined  systems  displayed  the  most  uniform  pressure  drop 
fluctuations  for  all  particle  sizes.  In  all  cases,  there  was  an  increase  in  pressure  drop 
fluctuation  with  decrease  in  gas  velocity  to  a  maximum.  Further  decrease  in  gas 
velocity  caused  a  decrease  in  fluctuation.  The  only  variation  to  this  behavior  was  for 
the  125pm  glass  beads  at  the  lower  solids  mass  flow  rate.  In  this  case  the  pressure 
drop  fluctuations  decreased  to  a  minimum,  then  increased  to  a  maximum,  and  then 
decreased  again  as  the  gas  velocity  was  decreased.  The  results  for  the  inclined  systems 
are  shown  in  Figures  A-34.  A-35  and  A-36. 

4.4.4  Electrostatic  Pressure  Drop 

The  effect  of  electrostatic  forces  on  the  pressure  drop  was  not  as  expected.  An 
increase  in  pressure  drop  was  anticipated  for  all  vertical  systems;  however,  some  cases 
resulted  in  pressure  drop  reduction.  The  horizontal  systems  also  showed  pressure  drop 
reductions  for  some  cases.  Only  the  45°  inclined  systems  showed  consistent  pressure 
drop  increase  with  electrostatics.  An  explanation  for  this  behavior  lies  in  the  particle 
flow  conditions.  Ally10  considered  the  particles  in  his  system  to  occupy,  on  an  average, 
an  annular  region  centered  at  D  /4  and  calculated  the  work  required  to  move  them  to 
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the  wall.  This  assumption  is  valid  for  a  homogenous  distribution  of  particles  moving 
with  negligible  radial  velocity  components.  This  condition  is  approached  when  D  is 
small  Observations  during  the  experiments  revealed,  for  some  cases,  a  considerable 
amount  of  internal  motion  in  the  flow.  This  apparently  changes  the  work  function  due 
to  electrostatics  by  altering  the  magnitude  or  direction  of  the  eletrical  field. 


A  comparison  was  made  between  the  inclined  systems  which  showed  predominantly 
positive  increases  to  the  pressure  drop  with  the  results  of  Ally30.  Ally  found  that 
electrostatic  pressure  drop  increased  with  decreasing  water  to  solids  ratio.  A  minimum 
water  to  solids  ratio  of  about  0.1  was  determined  above  which  there  was  no  appreciable 
increase  in  pressure  drop  due  to  electrostatics.  Figures  A-37,  A-38  and  A-39  show  the 
electrostatic  pressure  drop  vs.  water  to  solids  mass  ratio  for  the  inclined  systems.  For 
these  systems  the  minimum  water  to  solids  mass  ratio  is  between  0.03  and  0.07  which 
agrees  within  10%  of  Ally’s  result.  The  decrease  in  the  minimum  is  probably  due  to 
the  cleanliness  of  the  particles;  Ally's  being  cleaner. 

4.5  STABILITY 

The  stability  of  the  systems  studied  varied  widely  with  particle  size,  pipe  orien¬ 
tation  and  electrostatics.  In  the  vertical  systems,  motions  in  the  radial  and  tangential 
directions  in  the  pipe  became  very  pronounced  as  the  system  was  brought  to  choking. 
The  horizontal  systems  went  through  a  separation  across  the  radius  of  the  pipe,  with 
the  bottom  of  the  pipe  having  a  higher  solids  concentration.  In  the  45°  inclined  pipe, 
the  formation  of  retrograde  dunes  formed  for  all  cases. 

4.5.1  Choking 

The  phenomenon  of  choking  has  already  been  described  as  the  result  of  a  range 
of  instabilities.  It  is  probably  better  described  by  a  range  of  parameters.  For  this 
study,  the  choking  point  was  chosen  as  the  point  where  solids  first  dropped  out  of  the 
flow  below  the  feed  point  and  were  not  picked  up  by  the  gas  stream.  The  instabilities 
leading  to  this  point  occurred  before  this,  and  small  changes  in  the  system  operation 
could  have  brought  about  choking  at  different  times. 

Two  different  correlations  were  used  to  predict  the  choking  point.  They  were  the 
Yang'  and  the  Rose  and  Duckworth*6.  The  comparison  of  these  two  correlations  with 
the  experimental  results  is  shown  in  Table  4-2.  Both  correlations  ate  seen  to  under- 


predict  the  experimental  choking  point.  The  Rose  and  Duckworth  correlation  has  an 
error  range  of  831%  to  147%,  which  is  clearly  unsuitable  for  predicting  choking.  The 
Yang  correlation  does  much  better  with  errors  from  398%  to  26%.  The  largest  errors 
occur  for  the  79ym  glass  beads  which  apparently  act  as  clusters. 


The  effect  of  electrostatics  on  choking  varies  with  particle  size  and  mass  flow 
rate.  For  the  79ym  glass  beads,  electrostatics  tend  to  decrease  the  choking  gas  velocm 
for  the  lower  mass  flow  rates.  This  could  be  due  to  the  breaking  up  of  clusters  in  the 
system  due  to  the  like  charges  on  the  particles.  For  the  125ym  glass  beads  and  the 
128  ym  Plexiglas  beads,  the  presence  of  electrostatics  increases  the  choking  gas  velocity. 
The  effects  of  electrostatics  on  the  450ym  glass  beads  is  only  a  slight  increase  in  the 
choking  gas  velocity. 


4.5.2  Saltation 


Saltation  in  the  horizontal  pipe  was  determined  at  the  superficial  gas  velocity 
where  solid  particles  were  in  constant  contact  with  the  bottom  of  the  pipe.  The  450ym 
glass  beads  did  not  salt  out  in  this  study  as  the  vertical  section  leading  to  the  horizon¬ 
tal  section  became  plugged  before  any  saltation  occurred. 


Two  correlations  w'ere  compared  to  the  experimental  results.  The  correlation  of 
Owens*'  which  relates  a  pseudc  Froude  number  to  saltation  as  follows: 


P  17*  f  /2p  gD  >  0.01 

1  P*  (  r  r 


The  second  correlation,  is  that  of  Rizk‘  which  is 
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Both  of  these  expressions  are  implicit  in  saltation  gas  velocity  and  can  be  solved  by 
iteration.  The  results  of  these  calculations  and  the  experimental  saltation  velocities  are 
shown  in  Table  4-3. 

Table  4-3:  Saltation  Velocity  Analysis 

Part  leal  W  («xp«r)  u  ,  (axpar)  U  ,  (calc) 

a  gsalt  ^  gsalt 

-3 

yin  (kg/sxlO  High  R.H.  Low  R  H.  Ova  ns  Rlzk 


125 

18  3 

5  78 

4.36 

2.60 

5  24 

glass 

beads 

27.0 

5  92 

4  88 

2.60 

5  83 

79 

17.4 

5  38 

8  05 

2  00 

5  09 

glass 

beads 

26  5 

5.83 

6.  17 

2  00 

5.73 

450 

20  6 

• 

• 

5  41 

5  89 

glass 

baads 

30  5 

• 

• 

5  41 

6.50 

NO  SALTATION  OCCURRED 

For  the  79ym  glass  beads  the  Rizk  correlation  was  far  superior  to  the  Owens  correla¬ 
tion  for  both  high  and  low  humidity,  with  errors  of  2%  to  6$  for  the  high  humidity 
and  to  58%  for  the  low  humidity  The  Owens  correlation  was  off  by  303%  to 


16^% 


For  the  125 um  glass  beads  the  Rizk  correlation  was  again  superior.  Errors  ranged 
from  2%  to  17%,  while  the  Owens  correlation  ranged  from  68%  to  128%.  The  Owens 
correlation  was  better  for  the  125ym  particles  than  for  the  79ym  particles  but  still  was 
less  accurate  than  the  Rizk  correlation. 

The  450ym  glass  beads  had  no  experimental  saltation  velocities  as  the  particles 
never  salted  out.  Comparison  of  the  saltation  velocities  obtained  by  the  two  correla¬ 
tions  with  the  experimental  choking  velocities  for  450ym  glass  beads  shown  in  Table 

4-2  shows  that  the  expected  saltation  velocities  to  be  below  that  of  choking  which  is 

what  was  found. 

4.5.3  Linear  Analysis 

A  linear  stability  analysis  was  performed  for  all  conditions.  Values  for  the  solids 
friction  factor  were  obtained  by  solving  equation  2-15  for  f  using  the  experimental 
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pressure  drop,  gas  velocity,  and  particle  velocity.  Equation  2-31  was  then  solved  for 

m  and  m  .  For  all  cases,  m  was  negative  and  much  greater  than  m  .  The  condition 

12:  i 

for  stability  was  then  the  sign  of  m  For  positive  values  of  m  ,  the  flow  was  con- 

i  i 

sideied  unstable.  As  the  effect  of  the  electrostatic  forces  did  not  always  fit  the  form 
of  equation  2-25,  the  low  humidity  cases  were  treated  identically  to  the  high  humidity 
cases,  with  the  electrostatic  effects  being  combined  in  the  solids  friction  factor. 

a.  Vertical  Systems.  The  results  of  the  linear  analysis  for  the  79ym  glass  beads  in  the 
vertical  pipe  are  included  in  tables  B-20  through  B-25.  The  analysis  shows  this  particle 
to  be  very  unstable,  except  at  the  higher  gas  velocities.  There  is  considerable  fluctua¬ 
tion  in  the  eigenvalue  m  This  fluctuation  was  observed  in  the  pulsing  and  twisting 
motions  of  the  flow,  although  complete  plugging  of  the  pipe  did  not  occur  at  the  un¬ 
stable  gas  velocities  resulting  from  the  analysis.  The  effects  of  electrostatics  was  found 
to  be  minor,  except  in  the  highest  mass  flow  rate  condition.  In  this  case,  electrostatics 
caused  instability.  at  higher  gas  velocities. 

The  results  for  the  125yrr.  glass  beads  are  included  in  tables  B-6  through  B-I3 


This  particle  size  showed  unusual  behavior  with  change  in  mass  flow  rate  and  electros¬ 
tatics.  Stability  increased  with  increasing  mass  flow  rate.  At  low  mass  flow  rates, 
electrostatics  increased  stability.  As  mass  flow  rate  increased,  this  effect  was  damped 


The  results  for  the  450ym  glass  beads  are  included  in  tables  B-14  through  B-19. 
A  very  unusual  result  was  obtained  for  the  lowest  mass  flow  rate.  In  that  case, 
stability  was  seen  to  increase  from  unstable  to  stable  with  a  decrease  in  the  gas 
velocity. 

The  results  for  the  128ym  Plexiglas  beads  are  included  in  tables  B-26  through 

B-29.  The  128ym  Plexiglas  beads  were  seen  to  show  unusually  unstable  behavior  by 
linear  analysis.  Stability  was  predicted  for  the  higher  gas  flow  rates,  but  unstable 
predictions  occured  for  most  of  the  lower  gas  flow  rates.  Again,  observations  of  the 
Plexiglas  beads  during  flow,  pointed  out  the  twisting  and  pulsing  instabilities  without 
plugging  of  the  pipe. 

b.  Horizontal  Systems.  The  linear  analysis  results  for  ail  horizontal  systems  are  in¬ 
cluded  in  tables  B-34  through  B-41.  Stable  eigenvalues  were  obtained  for  all  con¬ 
ditions.  This  result  does  not  reflect  the  true  nature  of  the  flow  conditions.  The 

voidage  used  in  calculations  was  for  a  uniformly  dispersion  of  particles  across  the  pipe. 
This  was  not  the  true  condition  as  gas  flow  rate  was  decreased  as  the  flow  separated  to 
cause  an  increased  concentration  in  the  bottom  of  the  pipe.  The  equations  for  the 
analysis  are  incapable  of  predicting  this  gravity  effect.  A  possible  solution  would  be  to 
use  equation  2-16.  This  possibility,  however,  requires  knowledge  of  the  individual  fric¬ 
tion  factors,  which  were  not  obtained. 

c  Inclined  Systems  The  linear  analysis  results  for  the  inclined  systems  are  included 

in  tables  B-50  through  B-49.  The  results  show  an  increase  in  stability  with  decrease  in 
humidity  for  all  systems.  The  19ym  glass  beads  had  unstable  eigenvalues  for  the  high 
humidity  cases  only,  while  the  125ym  and  450ym  glass  beads  had  stable  eigenvalues  for 
a  cases  In  general,  stability  was  increased  with  increase  in  gas  velocity  and  increase 
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5.0  CONCLUSIONS 

1.  The  use  of  electrostatic  ring  probes  was  found  to  work  reasonably  well  in  the 
0.0508  m  pipe.  Agreement  with  correlations  tested  was  within  50%.  The  use 
of  this  type  probe  in  the  horizontal  and  inclined  orientations  is  not  recom¬ 
mended  in  the  probes  current  form  as  a  particle  density  gradient  can  exist  in 
these  orientations  which  can  cause  significant  deviations  in  the  signals  from 
the  probes. 

2.  Various  flow  patterns  were  observed  for  the  different  orientations  studied: 

a.  In  the  vertical  orientation,  radial  and  tangential  disturbances  were  ob¬ 
served  as  the  system  was  brought  toward  choking. 

b.  In  the  horizontal  systems,  the  particle  density  in  the  lower  portion  of 
the  pipe  increased  as  saltation  was  approached.  Prior  to  uniform  sal¬ 
tation,  blunt  nosed  dunes  formed  on  the  bottom  of  the  pipe. 

c.  In  the  inclined  orientation,  retrograde  dunes  formed  and  were  observed 
to  flow  along  the  bottom  of  the  pipe  against  the  main  stream. 

3.  The  pressure  drop  vs.  gas  velocity  was  found  to  depend  heavily  on  pipe 


orientation  and  panicle  size.  The  vertical  pipe  was  found  to  control  the  sys¬ 
tem  in  terms  of  stability  for  the  largest  particles.  The  horizontal  section  was 
found  to  control  the  system  for  the  smaller  panicles.  For  the  medium  par¬ 
ticle  size,  both  the  vertical  and  horizontal  sections  contributed  to  instabilities. 

4.  A  unique  behavior  w.as  found  for  the  125  pm  glass  beads.  The  pressure  drop 
curve  passed  through  a  maximum  as  gas  velocity  was  decreased.  This  be¬ 
havior  was  compared  to  similar  observations  made  by  Zaltash29who  attributed 
this  phenomenon  to  a  particle  to  tube  size  ratio.  As  this  behavior  was  not 
found  for  the  128ym  Plexiglas  particles,  a  density  factor  is  believed  to  con¬ 
tribute.  For  Zaltash's  study,  the  group  (D  /D)(p  /  p )  equaled  6.06.  For  this 

p  i  p  r 

study,  this  group  equaled  4.79. 

5.  Fluctuations  in  the  pressure  drop  were  found  to  decrease  with  increasing 
solids  mass  flow  rate.  For  the  smallest  particles,  there  was  little  change  in 
the  fluctuations  with  change  in  gas  velocity.  As  particle  size  increased,  the 
dependency  of  the  fluctuations  on  gas  velocity  increased,  with  fluctuations  in¬ 
creasing  with  decreasing  gas  velocity.  The  effect  of  electrostatics  was  found 
to  vary  considerably  with  particle  size  and  pipe  orientation. 

6.  The  effect  of  electrostatic  forces  on  the  systems  studied  was  found  to  vary 
considerably.  The  orientation,  particle  size,  and  gas  velocity  were  found  to 
affect  the  impact  of  electrostatics  significantly. 


7.  The  correlation  of  Konno  and  Saito  was  found  to  be  comparable  to  that  of 
Yang  in  the  vertical  systems  studied  with  the  advantage  of  over-predicting 
for  most  conditions.  In  the  horizontal  systems,  the  correlation  of  Yang  for 
this  orientation  was  found  superior. 

8.  The  two  correlations  used  to  predict  choking  gas  velocities  were  found  to  be 
inadequate  as  the  velocities  were  below  the  experimental  velocities  by  as  much 
as  831%.  The  use  of  the  Yang  correlation  for  choking  velocity  was  found  to 
be  superior  to  that  of  Rose  and  Duckworth,  however,  both  correlations 
under-predicted  the  choking  point  The  instabilities  associated  with  choking 
cover  a  wide  range  of  gas  velocities,  and  therefore,  the  definition  of  choking 
should  include  the  capacity  to  cover  this  range. 

9.  For  predicting  saltation,  the  correlation  of  Rizk  was  found  to  predict  within 
17%  of  the  experimental  values. 

10.  The  linear  analysis  performed  was  found  to  describe  the  vertical  systems  as 
far  as  observations  made.  Its  applicability  to  the  horizontal  and  inclined 
orientations  is  questionable  as  the  fundamental  equations  used  do  not  ade¬ 


quately  describe  the  system. 


6.0  RECOMMENDATIONS 

The  following  recommendations  are  included  as  possibilities  for  further  investiga¬ 
tion: 

1.  A  means  of  verifying  the  velocities  obtained  with  the  electrostatic  probes 
should  be  examined.  This  could  be  by  such  means  as  high  speed  photog¬ 
raphy. 

2.  A  selection  of  particles  with  different  densities  and  sizes  should  be  studied. 

In  conjunction  with  this,  other  pipe  sizes  should  be  included.  A  particle  size 

and  density  range  should  include  values  of  the  group  (D/D)(p  /p)  within 

r  '  p  f 

the  range  of  the  unique  behavior  found. 

3.  The  role  of  electrostatics  requires  much  further  attention.  Work  in  this  area 
should  attempt  to  evaluate  the  direction  and  magnitude  of  the  electric  field. 
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Figure  A-l:  Particle  Velocity  vs.  Gas  Velocity  for 

79vim  Glass  Beads  in  the  Vertical  Orientation 
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Figure  A -2:  Particle  Velocity  vs.  Gas  Velocity  for 
450iim  Glass  Beads  in  the  Vertical 
Orientation 


AV  AAlAA  A A.A  AAwj _  »  -A  «N  -A  »V  *- * . 


OCITY  (M/SEC) 


< 


Figure  A-4:  Partical  Velocity  vs.  Gas  Velocity  for 
79ym  Glass  Beads  in  the  Horizontal 
Orientation 
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Figure  A -5:  Partical  Velocity  vs.  Gas  Velocity  for 
450vitn  Glass  beads  in  the  Horizj  ital 
Orientation 
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Pressure  Drop  vs.  Gas  Velocity  for  the 
450ym  Glass  Beads  in  the  Inclined 
Orientation 
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Figure  A -16:  Pressure  Drop  vs.  Gas  Velocity  for 

79ym  Glass  Beads  for  the  Combined 
Orientations  at  tbe  Lower  Mass  Flow  Rate 
With  Electrostatics 
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Figure  A-17:  Pressure  Drop  vs.  Gas  Velocity  for 
79ym  Glass  Beads  for  the  Combined 
Orientations  at  the  Higher  Mass  Flow  Rate 
Without  Electrostatics 
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Figure  A-18:  Pressure  Drop  vs.  Gas  Velocity  for 
79ym  Glass  Beads  for  the  Combined 
Orientations  at  the  Higher  Mass  Flow  Rate 
With  Electrostatics 
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Figure  A -20:  Pressure  Drop  vs.  Gas  Velocity  for 
125ym  Glass  Beads  for  the  Combined 
Orientations  at  the  Lower  Mass  Flow  Rate 
With  Electrostatics 
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Figure  A -21:  Pressure  Drop  vs.  Gas  Velocity  for 
125ym  Glass  Beads  for  the  Combined 
Orientations  at  the  Higher  Mass  Flow  Rate 
Without  Electrostatics 
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Figure  A -25:  Pressure  Drop  vs.  Gas  Velocity  for 
450ym  Glass  Beads  for  the  Combined 
Orientations  at  the  Higher  Mass  Flow  Rate 
Without  Electrostatics 
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Figure  A -29:  Pressure  Drop  Fluctuation  vs. 

Gas  Velocity  for  450um  Glass  Beads 
in  the  Vertical  Orientation 
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Figure  A -32;  Pressure  Drop  Fluctuation  vs. 

Gas  Velocity  for  79ym  Glass  Beads 
in  the  Horizontal  Orientation 
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Figure  A -38:  Electrostatic  Pressure  Drop 

vs.  kg  H  O/kg  solids  for  the 

125ym  Glass  Beads  in  the  Inclined 
Orientation 
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Ratio  of  Experimental  to  Calculated 
Pressure  Drop  vs.  Gas  Velocity  for 
79ym  Glass  Beads  in  the  Vertical 

Orientation  with  W  =18.  xlO  3 
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Figure  A-42;  Ratio  of  Experimental  to  Calculated 
Pressure  Drop  vs.  Gas  Velocity  for 
79ym  Glass  Beads  io  the  Vertical 

Orientation  with  W  =27.  xlO  3 
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Figure  A -43:  Ratio  of  Experimental  to  Calculated 
Pressure  Drop  vs.  Gas  Velocity  for 
125vm  Glass  Beads  in  the  Vertical 

Orientation  with  W  =8.  xlO'3 
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Figure  A-45:  Ratio  of  Experimental  to  Calculated 
Pressure  Drop  vs.  Gas  Velocity  for 
125pm  Glass  Beads  in  the  Vertical 

Orientation  with  W  =27.  xlO  3 
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Figure  A -50:  Ratio  of  Experimental  to  Calculated 

Pressure  Drop  vs.  Gas  Velocity  for 
128pm  Plexiglas  Beads  in  the  Vertical 

Orientation  with  W  =12.  xlO  3 
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Figure  A-51:  Ratio  of  Experimental  to  Calculated 
Pressure  Drop  vs.  Gas  Velocity  for 
79ym  Glass  Beads  in  the  Horizontal 

Orientation  with  W  si 7.  xlO’3 
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Figure  A -52: 

Ratio  of  Experimental  to  Calculated 
Pressure  Drop  vs.  Gas  Velocity  for 

79  vim  Glass  Beads  in  the  Horizontal 
Orientation  with  W  =24.  xlO  3 
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Figure  A-54:  Ratio  of  Experimental  to  Calculated 
Pressure  Drop  vs.  Gas  Velocity  for 
125ym  Glass  Beads  in  the  Horizontal 

Orientation  with  W  =24.  xlO  3 


Table  B-2:  Comparison  of  Absolute  Mean  Percent 

Error  between  the  Correlation  of  Konno  and 
Saito  and  that  of  Yang  for  125iim  Glass 
Beads  in  the  Vertical  Orientation 
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Table  B-3:  Comparison  of  Absolute  Mean  Percent 

Error  between  the  Correlation  of  Konno  and 
Saito  and  that  of  Yang  for  450ym  Glass 
Beads  in  the  Vertical  Orientation 


w 

s 

Humidity 

Konno  and 

Saito 

Yang 

(x  10 

3) 

Abs.  Mean 

Sid 

Abs.  Mean  Std 

Eiror  (%) 

Dev 

Error 

(%)  Dev 

12.5 

High 

21.7 

12.1 

36.5 

11.7 

10.0 

Low 

25.7 

13.9 

34.8 

16.5 

22.7 

High 

12.2 

6.9 

37.1 

10.5 

20.6 

Low 

32.3 

22.0 

23.9 

14.6 

33.3 

High 

11.5 

10.3 

24.2 

14.1 

30.0 

Low 

19.3 

7.4 

16.5 

11.6 

Table  B-4:  Comparison  of  Absolute  Mean  Percent 

Error  between  the  Correlation  of  Konno  and 
Saito  and  that  of  Yang  for  128 uni  Plexiglas 
Beads  in  the  Vertical  Orientation 

Humidity  Konno  and  Saito  Yang 

Abs.  Mean  Std  Abs.  Mean  Std 

Error  (%)  Dev  Error  (%)  Dev 


Table  B-5:  Comparison  of  Absolute  Mean  Percent 

Error  between  the  Correlation  of  Konno  and 
Saito  and  that  of  Yang  for  79ym  Glass 
Beads  in  the  Horizontal  Orientation 


W  Humidity  Konno  and  Saito  Yang 


(X  10") 

Abs.  Mean 

Error  (%) 

Std 

Dev 

Abs.  Mean  Std 

Error  (%)  Dev 

19.0 

High 

17.6 

11.0 

12.4 

12.9 

15.4 

Low 

10.2 

7.5 

13.2 

7.9 

24.1 

High 

12.1 

3.9 

24.1 

10.6 

24.5 

Low 

14.8 

5.1 

14.0 

6.9 

Table  B-6:  Comparison  of  Absolute  Mean  Percent 

Error  between  the  Correlation  of  Konno  and 
Saito  and  that  of  Yang  for  125}ini  Glass 
Beads  in  the  Horizontcal  Orientation 


W  Humidity  Konno  and  Saito  Yang 


(x  10'3) 


Abs.  Mean  Std  Abs.  Mean  Std 
Error  (%)  Dev  Error  (%)  Dev 


18.5 

High 

29.2 

9.3 

10.1 

8.6 

15.8 

Low 

19.8 

12.9 

24.1 

16.4 

25.3 

High 

55.9 

3.8 

36.7 

2.8 

23.5 

Low 

29.2 

11.5 

9.9 

8.9 

Table  B-7:  Comparison  of  Absolute  Mean  Percent 

Error  between  the  Correlation  of  Konno  and 
Saito  and  that  of  Yang  for  450ym  Glass 
Beads  in  the  Horizontal  Orientation 

Humidity  Konno  and  Saito  Yang 


Abs.  Mean  Std  Abs.  Mean  Std 

Error  (%)  Dev  Error  (%)  Dev 
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Table  B-ll:  Eigenvalues  for  125^01  Glass  Beads 
in  the  Vertical  Orientation 

When  Ws=l8.7xl0  3  and  ILH.=  52.0 


RUN  U  U  m  m 

g  f  l  2 


19 

16.76 

12.23 

-2.725 

-64.098 

20 

13.38 

7.59 

-1.385 

-92.192 

21 

10.63 

10.16 

-1.747 

-30.931 

22 

8.36 

6.17 

-0.885 

-71.308 

23 

6.14 

5.04 

-1.849 

-61.962 

24 

5.80 

5.74 

-1.426 

-19.581 

25 

5.02 

4.68 

-1.278 

-41.258 

76 

10.40 

8.36 

-1.577 

-58.488 

77 

8.22 

7.86 

-1.728 

-31.751 

78 

7.45 

5.69 

-0.692 

-69.590 

79 

6.78 

6.74 

-1.411 

-16.090 

80 

6.49 

4.97 

-0.632 

-71.135 

81 

5.78 

5.69 

-1.438 

-22.182 

iU 


Table  B-13:  Eigenvalues  for  125yin  Glass  Beads 
in  the  Vertical  Orientation 


Table  B-14:  Eigenvalues  for  450vim  Glass  Beads 
in  the  Vertical  Orientation 

When  Ws  =  12.5xl0  3  and  R-H.=  50.0 


9.71 

1.376 

-7.554 

9.71 

2.500 

-6.998 

5.62 

1.472 

-9.199 

6.41 

0.099 

-6.172 

4.89 

-0.693 

-8.255 

5.04 

-1.643 

-7.182 

4.49 

-2.326 

-8.181 

8.36 

1.001 

-7.321 

5.55 

0.083 

-8.566 

Table  B-17:  Eigenvalues  for  450um  Glass  Beads 
in  the  Vertical  Orientation 

When  Ws  =  20.6xl0  3  and  R-H.=  16.0 


RUN 

u 

s 

u 

f 

m 

i 

i 

129 

15.34 

12.01 

4.012 

-10.039 

130 

14.00 

11.79 

-1.009 

-9.952 

131 

12.36 

11.19 

-0.674 

-8.801 

132 

10.74 

7.26 

-0.281 

-13.441 

133 

10.05 

8.36 

-0.068 

-10.296 

134 

9.49 

5.55 

0.754 

-15.969 

135 

9.16 

4.37 

2.508 

-19.249 

136 

7.87 

3.86 

1.173 

-19.938 

137 

6.58 

3.51 

1.059 

-18.505 

138 

5.89 

2.75 

1.562 

-19.867 

139 

16.58 

14.05 

-1.462 

-9.934 

Table  B-20:  Eigenvalues  for  79ym  Glass  Beads 
in  the  Vertical  Orientation 


When  Ws= 9.5x1 0‘3  and  R-H.=  53.7 

RUN 

*  u 

& 

u 

f 

m 

i 

m 

2 

150 

11.16 

8.36 

-0.820 

-100.078 

151 

9.11 

6.11 

-0.087 

-118.126 

152 

8.36 

3.86 

1.164 

-172.593 

153 

7.62 

3.51 

3.612 

-175.384 

154 

6.62 

3.77 

2.181 

-146.208 

155 

5.82 

2.29 

4.985 

-213.789 

156 

6.53 

5.69 

0.731 

-73.453 

157 

5.20 

3.07 

4.568 

-143.880 

158 

4.18 

3.37 

2.986 

-80.709 

159 

4.71 

3.93 

1.686 

-81.802 

Table  B-21:  Eigenvalues  for  79um  Glass  Beads 
in  the  Vertical  Orientation 

When  Ws=8.8xl0  3  and  R-H.=  16.1 


RUN 

u 

i 

u 

f 

m 

i 

r 

200 

12.00 

5.00 

0.226 

-171.283 

201 

10.43 

6.60 

-0.243 

-119.470 

202 

8.31 

6.60 

-0.319 

-88.570 

203 

7.16 

4.89 

1.171 

-112.449 

204 

5.69 

5.16 

1.524 

-60.401 

205 

5.91 

4.23 

2.110 

-107.413 

206 

5.38 

3.83 

2.232 

-108.191 

207 

5.10 

2.71 

4.700 

-155.089 

208 

4.36 

3.83 

2.422 

-60.254 

209 

6.94 

5.74 

-0.718 

-82.480 

/.v.v.v 
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Table  B-23:  Eigenvalues  for  79vtm  Glass  Beads 
in  the  Vertical  Orientation 

When  Ws= 18.3x1  O'3  and  R-H.=  16.0 

RUN 

u 

s 

u 

f 

m 

i 

m 

190 

14.34 

12.01 

-2.105 

-102.406 

191 

12.36 

6.67 

-0.738 

-192.559 

192 

10.85 

5.24 

1.031 

-221.892 

193 

10.67 

5.04 

0.335 

-237.665 

194 

8.31 

6.17 

-1.103 

-103.003 

195 

9.20 

6.41 

0.843 

-144.561 

1% 

7.54 

5.37 

0.082 

-141.464 

197 

6.47 

6.11 

-0.138 

-69.074 

198 

5.56 

5.00 

1.073 

-89.761 

.  V.  •  •  ■  .  •  .  •  ■ 
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Table  B-27:  Eigenvalues  for  128um  PLEXIGLAS  BEADS 

in  the  Vertical  Orientation 

When  Ws=  8.8x10  3  and  R-H.=  16.0 


u 

t 

u 

f 

m 

i 

t 

9.82 

6.17 

-1.300 

-104.704 

8.45 

5.50 

-0.117 

-98.514 

6.14 

4.75 

0.632 

-80.790 

4.62 

3.25 

2.279 

-94.641 

3.55 

2.43 

2.248 

-92.672 

4.27 

2.90 

3.012 

-102.794 

3.32 

2.95 

-0.034 

-54.146 

8.34 

7.18 

-1.064 

-60.789 

6.53 

5.50 

1.044 

-64.348 

.  .  *  4  -  .  *r-  '  ./  ■ 
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Table  B-28:  Eigenvalues  for  128um  PLEXIGLAS  BEADS 

in  the  Vertical  Orientation 

When  Ws=12.9xl0  3  and  R.H.=  52.8 


RUN  U  U  m  m 

t  f  i  : 


240 

12.97 

6.88 

-2.667 

-137.974 

241 

12.63 

7.96 

-1.694 

-121.708 

242 

8.22 

7.59 

-1.375 

-56.763 

243 

6.00 

3.92 

1.055 

-126.422 

244 

4.74 

3.12 

-0.704 

-127.767 

245 

8.15 

5.20 

0.263 

-126.677 

246 

3.63 

3.18 

-0.620 

-76.412 

248 

8.04 

5.95 

-0.897 

-105.060 

249 

4.89 

3.18 

-0.288 

-136.043 

jflrfi rfvY . 
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Table  B-30:  Eigenvalues  for  125um  Glass  Beads 
in  the  Horizontal  Orientation 


When  Ws= 

18.5xl0'3  and  R-H.=  51.7 

u 

u 

m 

m 

1 

f 

i 

2 

i 

7,26 

-4.575 

-80.759 

! 

6.67 

-3.805 

-87.783 

i 

8.69 

-5.980 

-51.978 

) 

6.11 

-4.665 

-87.534 

) 

3.24 

-4.448 

-160.639 

5 

3.53 

-3.920 

-133.033 

5 

3.46 

-3.644 

-123.811 

1 

2.33 

-3.188 

-164.382 

Table  B-30:  Eigenvalues  for  125ym  Glass  Beads 
in  the  Horizontal  Orientation 

When  Ws  =  18.5xl0  3  and  ILH.=  51.7 
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Table  B-31:  Eigenvalues  for  125vim  Glass  Beads 
in  the  Horizontal  Orientation 

When  Ws=15.8xl0  3  and  R.H.=  16.9 


RUN 

u 

s 

u 

r 

m 

i 

i 

255 

9.45 

3.93 

-1.614 

-129.698 

256 

9.20 

4.89 

-3.139 

-101.078 

257 

7.56 

3.24 

-3.103 

-135.064 

258 

6.49 

4.86 

-2.682 

-71.661 

259 

6.00 

2.90 

-2.578 

-124.840 

260 

4.36 

2.43 

-4.621 

-63.201 

261 

10.31 

7.59 

-3.718 

-70.016 

tiSt 
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Table  B-34:  Eigenvalues  for  79ym  Glass  Beads 
in  the  Horizontal  Orientation 

When  Ws=19.0xl0  3  and  R.H.=  55.4 


U  m 

g  f  i 


5.00 

-2.756 

-217.775 

7.96 

-3.980 

-116.339 

5.74 

-4.528 

-172.826 

5.69 

-3.945 

-171.717 

4.75 

-3.600 

-185.832 

4.13 

-2.536 

-205.899 
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Table  B-35:  Eigenvalues  for  79ym  Glass  Beads 
in  the  Horizontal  Orientation 

When  Ws=15.4xl0  3  and  R.H.=  19.1 


RUN 

u 

g 

u 

f 

m 

i 

r 

284 

11.60 

4.59 

-3.510 

-252.039 

285 

11.29 

3.37 

-3.454 

-327.862 

286 

10.63 

8.36 

-3.755 

-121.101 

287 

9.87 

7.03 

-3.277 

-144.643 

288 

9.20 

8.36 

-3.410 

-83.293 

289 

8.00 

4.49 

-2.711 

-204.325 

290 

5.87 

4.49 

-3.526 

-97.430 

_ 
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Table  B-36:  Eigenvalues  for  79ym  Glass  Beads 
in  the  Horizontal  Orientation 

When  Ws=24.1xl0  3  and  R.H.=  53.4 


5.00 

-3.459 

-267.717 

5.69 

-2.356 

-200.0% 

2.29 

-2.779 

-544.860 

2.75 

-2.520 

-434.440 

2.75 

-2.447 

-277.358 

7.18 

-2.813 

-112.160 

3.93 

-2.610 

-317.180 

Table  B-38: 

Eigenvalues  for  450ym  Glass  Beads 
in  the  Horizontal  Orientation 

When  Ws  =  19.8xl0  3  and  R.H.=  56.3 

RUN 

U  U 

t  f 

m  m 

1  2 

319 

14.09 

10.48 

-2.404 

-11.995 

320 

10.94 

7.26 

-3.569 

-14.038 

321 

10.36 

4.86 

-3.287 

-19.096 

322 

8.54 

5.90 

-3.443 

-13.878 

323 

7.69 

4.75 

-3.576 

-16.060 

324 

6.89 

3.44 

-3.031 

-20.199 

mjT^-LX  .tjCunV-* . J*kjf  ..-r  A~ jJjjUtULa.  ':?,^£L.-,\  t, l£A 


>*■ 


1  V*  ’•“  V  v  '.*  \fc  ’v  ■ 
-  *  ■*  * 


Table  B-43:  Eigenvalues  for  125ym  Glass  Beads 
in  the  Inclined  Orientation 

When  Ws  =  17.8xl0  3  and  R.H.=  19.0 


9.30 

-2.515 

-41.271 

5.37 

-1.999 

-102.630 

7.26 

-2.872 

-69.243 

5.95 

-2.061 

-74.621 

3.93 

-1.386 

-115.478 

4.89 

-3.531 

-74.407 

3.06 

-20.342 

-94.250 

-  *  -  -  ^ 


Table  B-44:  Eigenvalues  for  125ym  Glass  Beads 
in  the  Inclined  Orientation 

When  Ws=27.0xl0  3  and  R.H.=  56. 


Table  B-45:  Eigenvalues  for  125utn  Glass  Beads 
in  the  Inclined  Orientation 

When  Ws  =  28.9xl0  3  and  R.H.=  22.0 


^U)-fll84  074  GAS-SOL  ID  TRANSPORT  IN  A  00508  M  PIPE  AT  VARIOUS  2/2 

INCLINATIONS  UITH  AND  UITHOUT  ELECTROSTATICS(U)  ARHV 
HILITARV  PERSONNEL  CENTER  ALEXANDRIA  VA  C  A  NVLER 
UNCLASSIFIED  AUG  85  F/G  12/11  HL 
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Table  B-46:  Eigenvalues  for  450um  Glass  Beads 
in  the  Inclined  Orientation 

When  Ws=18.6xl0'3  and  R.H.=  55.7 

RUN 

U  U  m  m 

If  1  2 

347 

11.91 

9.86 

-2.581 

-10.295 

348 

11.11 

6.60 

-1100 

-14.630 

349 

10.82 

6.41 

-2.151 

-14.863 

350 

10.00 

6.17 

-1.925 

-14.417 

351 

9.29 

5.04 

-1.575 

-16.556 

352 

7.38 

4.72 

-2.457 

-15.117 

Table  B-47: 


Eigenvalues  for  450um  Glass  Beads 
in  the  Inclined  Orientation 


When  Ws=20.3xl0'3  and  R.H.=  21.9 


RUN 

U 

t 

U 

r 

m 

i 

r 

333 

10.89 

5.90 

-5.218 

-17.784 

334 

9.82 

4.59 

-3.273 

-19.685 

335 

9.29 

5.33 

-4.313 

-17.444 

336 

8.85 

5.37 

-3.969 

-16.101 

337 

8.16 

3.59 

-3.510 

-21.504 

339 

7.16 

4.37 

-6.960 

-19.263 

m 


Table  B-48: 

Eigenvalues  for  450ym  Glass  Beads 
in  the  Inclined  Orientation 

When  Ws=29.1xl0*3  and  R.H.=  53.6 

RUN 

U 

u 

m 

m 

t 

f 

l 

2 

354 

11.23 

6.67 

-2.868 

-18.353 

355 

10.61 

5.04 

-1.694 

-22.584 

356 

10.11 

7.96 

-2.450 

-13.285 

357 

8.72 

5.50 

-1.581 

-18.121 

358 

8.28 

4.03 

-1.379 

-24.286 

359 

7.72 

6.60 

-3.122 

-12.371 

172 


Table  B -49:  Eigenvalues  for  450ym  Glass  Beads 

in  the  Inclined  Orientation 

When  Ws=28.8xl0*3  and  R-H.=  18.0 


RUN  U  U  m  m 

g  f  1  2 
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Table  B-50:  Eigenvalues  for  79ym  Glass  Beads 
|  in  the  Inclined  Orientation 

When  Ws=18.1xl0'3  and  R.H.=  52.8 


RUN 

U 

t 

U 

r 

m 

i 

i 

403 

10.05 

7.18 

-2.216 

-148.197 

404 

9.78 

6.11 

-2.055 

-177.720 

405 

9.05 

5.20 

-2.100 

-197.538 

406 

8.45 

7.26 

-1.751 

-103.945 

407 

7.22 

2.71 

1.387 

-325.563 

408 

7.02 

2.39 

1.738 

-360.726 

409 

6.78 

3.35 

-0.454 

-251.050 
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Table  B-51:  Eigenvalues  for  79ym  Glass  Beads 
in  the  Inclined  Orientation 

When  Ws=16.5xl0"3  and  R.H.=  18.5 


RUN 

U 

s 

U 

r 

m 

i 

m 

389 

11.08 

8.81 

-2.401 

-99.369 

390 

10.36 

4.13 

-0.529 

-271.150 

391 

9.52 

5.69 

-1.205 

-183.498 

392 

8.05 

7.59 

-2.110 

-70.350 

393 

7.60 

6.11 

-1.646 

-123.537 

394 

6.80 

5.16 

-1.584 

-140.582 

175 


Table  B-52:  Eigenvalues  for  79nm  Glass  Beads 
in  the  Inclined  Orientation 

When  Ws= 26.6x1 0‘3  and  R.H.=  54.6 


RUN 

U 

i 

U 

r 

m 

i 

m 

2 

410 

9.11 

6.88 

-2.159 

-169.376  • 

411 

7.95 

5.95 

-0.662 

-179.888 

412 

7.06 

5.20 

-0.891 

-186.387 

413 

6.84 

3.44 

0.769 

-322.710 

415 

6.45 

3.12 

1.338 

-332.440 

416 

6.25 

4.75 

-0.817 

-179.629 

176 
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Table  B-53:  Eigenvalues  for  79ym  Glass  Beads 
in  the  Inclined  Orientation 


When  Ws= 25.8x1 0*3  and  R.H.=  17.6 


RUN 

U 

8 

U 

r 

m 

i 

i 

396 

9.22 

8.36 

-2.824 

-107.020 

397 

8.61 

6.17 

-2.313 

-194.539 

398 

8.78 

5.20 

-1.532 

-254.876 

399 

6.89 

4.89 

-1.717 

-200.346 

400 

6.06 

4.75 

-3.326 

-130.947 

401 

7.39 

4.75 

-2.912 

-239.469 

402 

7.64 

5.90 

-2.860 

-174.742 

m 


APPENDIX  C 
COMPUTER  PROGRAMS 


CR0S1.F0R 

CROSSCORRELATION  PROGRAM 

THIS  IS  THE  FIRST  PROGRAM  IN  A  SERIES  OF  TWO  WHICH 
MUST  BE  EXECUTED  IN  ORDER  TO  CALCULATE  THE  PARTICLE 
VELOCITY.  THIS  PROGRAM  TAKES  THE  ANALOG  SIGNAL  FROM 
THE  TAPE  RECORDER  AND  CONVERTS  IT  INTO  A  DIGITAL 
SIGNAL.  DATA  IS  INPUTED  FROM  THE  FIRST  TWO  CHANNELS 
ON  THE  A/D  BOARD.  THE  UPSTREAM  SIGNAL  SHOULD  BE 
CONNECTED  TO  CHANNEL  0  AND  THE  DOWNSTREAM  SIGNAL 
SHOULD  BE  CONNECTED  TO  CHANNEL  1  ON  THE  A/D  BOARD. 
DATA  SHOULD  BE  INPUTED  FOR  8  SECONDS  IN  ORDER  TO 
COLLECT  8000  POINTS  IN  FILE  ’FILE.DAT. 

THE  EXAMPLE  BEGINS  BY  DEFINING  CERTAIN  ARRAYS 
REQUIRED  BY  THE  SWEEP  PROCESS.  INFO  IS  USED  AS  THE 
SWEEP-INFORMATION  ARRAY  AND  IBUF(1,N+1)  IS  USED  AS 
THE  Nth  BUFFER. 

IN  ADDITION,  ARRAY  FBUF  IS  DEFINED.  THIS  ARRAY  IS 
USED  FOR  TEMPORARY  STORAGE  OF  A/D  CONVERTED  VALUES. 

DIMENSION  INFO(40),IBUF(1000,8),FBUF(1000) 

THE  FIRST  OPERATION  PERFORMED  IN  THIS  EXAMPLE  IS  TO 
CHECK  THE  STATUS  AND  READINGS  OF  A/D  CHANNELS  0  AND 
8.  NOTE  THAT  CHANNEL  0  CAN  NEVER  RETURN  ANY  STATUS 
OTHER  THAN  SINGLE-ENDED  VOLTAGE  INPUT  WITH  A  GAIN 
VALUE  OF  1.  CHANNEL  8,  ON  THE  OTHER  HAND.  MAY  HAVE  A 
PREAMPLIFIER  (THAT  IS,  AN  MNCAG)  ASSOCIATED  WITH  IT. 
THUS  ALLOWING  A  VARIABLE  STATUS  AS  WELL  AS  INPUT. 

THIS  SECTION  LOOPS  ON  OPERATOR  CONTROL  SO  THAT  * 
VARIATIONS  IN  INPUT  STATUS  AND  VALUE  CAN  BE 
VERIFIED. 

BEFORE  PROCEEDING,  SOME  ADDITIONAL  ARRAYS  NEED  TO  BE 
DEFINED. 

DIMENSION  RANGE(5,3),LENCK3),GAIN(5) 

INTEGER*4  IDOUBL 
REAL  *8  VALN 
LOGICAL*!  TYPE(10,3),ANS 


OUUU  L>  u  u  u  uuu 
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DATA  RANGE/.01..1.1.02.10.24.5.12,  1V0LTS 
.01..  1.1.02.10.24,0.0,  IMILLIAMPS 

.1.1.02,10.24,102.4.0.0/  !K  OHMS 

DATA  TYPE/*  7V’,’0’.’L7T7S’.4*’  ’, 

’  ’.’M’.T.’L’.’L’.T.’A’.’M’.’P’.’S’, 

'  7K7  707H7M7S’,3*’  7 
DATA  LENG/6,10,7/ 

DATA  GAIN /.5,5.,50.,500.,1. / 

DO  10  1=0, 8, 8 

NOW  GET  MNCAD  #0  STATUS  INFORMATION  FOR  THE 
APPROPRIATE  CHANNEL 

CALL  MADSTS(IND.I.INNER,IOUTER,NOMAG) 

NOW  READ  AND  CONVERT  THE  VALUE  RECEIVED  FROM  THE 
A/D  CHANNEL 

VAL=CAD2FP(IADINP(,I.IVAL)) 

NOW  NORMALIZE  THE  CONVERTED  VALUE  AND  SET  THE  DOUBLE 
PRECISION  RESULT  TO  A  FLOATING  POINT  NUMBER. 

CALL  KAD2DI(IVAL.IDOUBL) 

VALN=DJFLT(IDOUBL) 

NOW  PRINT  RAW  DATUM  (IV AL),  THE  NORMALIZED  DATUM 
(VALN).  AND  THE  ACTUAL  DATUM  (IV AL). 

TYPE  1000,I.IVAL,VAL,(TYPE(J,INNER),J=1,LENG(INNER)) 

.000  FORMAT!/'  THE  RAW  INPUT  READ  ON  CHANNEL  #’I1’  IS 

1  *06’  OCTAL.’/’  THE  ACTUAL  INPUT  IS  THEREFORE  ’F10.4, 

2  10A1) 

TYPE  1001,  VALN 

.001  FORMATC  THE  NORMALIZED  REPRESENTATION  OF  THE  INPUT 
1  IS  ’F10.07) 

REPORT  ON  THE  STATUS  OF  THE  CHANNEL 
TYPE  1004 

1004  FORMATC  THE  INPUT  IS  ’$) 

GO  TO  (14,15,16),NOMAG+2 
14  TYPE  1005 
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1005  FORMATO  QUASI-DIFFERENTIAL  ’$) 

GO  TO  30 

15  TYPE  1006 

1006  FORMAT!’  DIFFERENTIAL  ’$) 

GO  TO  30 

16  TYPE  1007 

1007  FORMAT!’  SINGLE-ENDED  ’$) 

30  GO  TO  (21.22.22.22, 22,22).IOUTER 

21  TYPE  1008 

1008  FORMAT!’  WITH  THE  GAIN  SET  TO  PROGRAMMABLE.’) 

GO  TO  10 

22  TYPE  1009,RANGE(IOUTER-1, INNER), (TYPE(J, INNER)  ,J-1, 

1  LENG!INNER)) 

1009  FORMAT!’  WITH  A  RANGE  OF  +  OR  -  ’F6.2.10A1) 

TYPE  1010,GAIN(IOUTER-1) 

1010  FORMAT!’  THE  EFFECTIVE  GAIN  IS  ’F5.1V) 

10  CONTINUE 

C 

C  CHECK  IF  THE  OPERATOR  WISHES  TO  VERIFY  INFORMATION 

C  AGAIN. 

C 

TYPE  1020 

1020  FORMAT!/’  SHOULD  CHANNELS  BE  CHECKED  AGAIN?  ’$) 
ACCEPT  2000,  ANS 
2000  FORMAT(Al) 

IF(ANS.EQ.'Y’)GO  TO  1 
C 

C  THE  EXAMPLE  NOW  GIVES  THE  OPERATOR  THE  OPTION  TO 

C  START  A  SWEEP  OF  CHANNEL  SWEEPS  ON  ANY  CHANNEL 

C  BETWEEN  0  AND  8.  UP  TO  A  MAXIMUM  OF  FOUR  CHANNELS 

C  MAY  BE  SELECTED  TO  BE  SAMPLED  DURING  EACH  CHANNEL 

C  SWEEP.  AN  AGGREGATE  TOTAL  OF  100  SAMPLES  WILL  BE 

C  TAKEN  IN  THIS  PROCESS. 

C 

C  THE  OPERATOR  MAY  NOW  ALSO  CHOOSE  THE  INTERSAMPLE 
C  INTERVAL  TIME,  OR  DWELL,  AS  A  INTEGRAL  MULTIPLE  OF 
C  .01  SECONDS.  THE  VALUE  CHOSEN  MUST  BE  LESS  THAN 

C  65536  BUT  GREATER  THAN  ZERO. 

C 

C  THE  OPERATOR  MAY  ALSO  CHOOSE  TO  HAVE  THE  INTERVAL  BE 

C  THE  TIME  BETWEEN  SUCCESSIVE  FIRINGS  OF  SCHMITT 

C  TRIGGER  1  OF  THE  PRIMARY  CLOCK  BY  CHOOSING  ZERO  (0) 

C  AS  THE  LENGHT  OF  THE  SAMPLE  PERIOD. 

C 


50  TYPE  1021 

1021  FORMATC  WHICH  CHANNEL  IS  THE  FIRST  A/D  CHANNEL  TO 

1  BE  SAMPLED  AT  THE  END  OF  EACH'/’  SAMPLE  PERIOD(0-8)? 

2  '$) 

ACCEPT  2001, ICHN 
2001  FORMATU5) 

IFUCHN.LT.O.OR.ICHN.GT.8)  GO  TO  50 

51  TYPE  1022 

1022  FORMATC  HOW  MANY  CHANNELS  SHOULD  BE  SAMPLED  EACH 
1  SAMPLE  PERIOD  (1-4)?  '$) 

ACCEPT  2001.NCHN 

IF(NCHN.LT.1.0R.NCHN.GT.4)  GO  TO  51 
C 

C  NOTE  THAT  THE  MODE  INDICATES  THAT  THE  SWEEP  WILL  NOT 

C  START  WHEN  THE  CALL  TO  THE  ADSWP  SUBROUTINE  IS  MADE 

C  IF  THE  SWEEP  IS  DRIVEN  BY  THE  PRIMARY  CLOCK  AT  A 

C  SPECIFIED  RATE.  THE  PROGRAM  WAITS  FOR  A  COMMAND  FROM 

C  THE  OPERATOR  AND  THEN  CALLS  THE  DIGO  SUBPROGRAM  TO 

C  START  THE  SWEEP  IMMEDIATELY. 

C 

C  IF  THE  SWEEP  IS  EXTERNALLY  DRIVEN,  i.e.,  BY  INPUT 

C  DIRECTLY  GIVE  TO  THE  MNCAD  MODULE.  OR  FROM  INPUT  TO 

C  ST2  OF  THE  PRIMARY  CLOCK  (MNCKW),  THE  PROGRAM  WAITS 

C  TILL  THE  USER  GIVES  THE  COMMAND  BEFORE  CALLING  THE 

C  SUBPROGRAM  ADSWP.  ADSWP  ARMS  AND  ENABLES  THE 

C  SAMPLING  PROCESS  IMMEDIATELY  IN  THIS  MODE. 

C 

TYPE  1019 

1019  FORMAT(5X,’SPECIFY  NUMBER  OF  CYCLES(1.2,3.4.5  OR  6) 

1  ’S) 

ACCEPT  2019.M 
2019  FORMAT(I) 

TYPE  2222 

2222  F ORM AT (5X, ’ WHI CH  UNIT  IS  TO  BE  OPENED?’  $) 

ACCEPT  2223.IUNIT 

2223  FORMAT(I) 

IF(IUNIT.EQ.3)  GO  TO  666 

OPEN(UNIT=2,NAME=’FILE.DAT\TYPE=’NEW\DISP=’SAVE’) 

GO  TO  667 

666  OPEN(UNIT=3,NAME=’FBUFl.DAT’,TYPE=’NEW’,DISP=’SAVE’ 

1  ) 

667  DO  31  I=1,M 

CALL  SETIBF(INFO,IND,IBUF(l,l),IBUF(1.2),IBUF(1.3), 

1  IBUF(1,4),IBUF(1,5),IBUF(1,6),IBUF(1,7),IBUF(1,8)) 


CALL  RLSBUF(INFO.IND.0.1.2,3.4.5,6,7) 

CALL  CL0CKA(4„IND) 

MODE=256 

CALL  ADSWP(INFO,100G,8,MODE.-1„.ICHN,NCHN) 

IF(INFO(1).NE.O)  STOP  ’ADSWP  -  LIBGEN  OR 
1  CONFIGURATION  ERROR.’ 

C 

C  THE  DEFAULTED  ARGUMENTS  INDICATE  THAT  THE  SWEEP  HAS 
C  NO  DELAY  FROM  THE  START  EVENT,  AND  NO  COMPLETION 

C  ROUTINE  IS  ASSOCIATED  WITH  IT. 

C 

C  AT  THE  OPERATOR’S  COMMAND  THE  SWEEP  WILL  BEGIN. 

C 

TYPE  1024 

1024  FORMATO  STRIKE  THE  "RETURN"  KEY  TO  START  THE  SWEEP 
1  .  ’$) 

ACCEPT  2000. ANS 
CALL  DIGOL..IND) 

C 

C  WAIT  FOR  EACH  BUFFER,  CONVERT  THE  INPUT  TO  ACTUAL 

C  VALUES,  TYPE  THE  ACTUAL  VALUES,  AND  THEN  GO  GET 

C  NEXT  BUFFER. 

C 

60  CALL  CAD2FP(IBUF(l.(IWTBUF(INFO..ID.INDM)).FBUF. 

1  1000) 

IF(IUNIT.EQ.3)  GO  TO  444 
WRJTE(2,33)FBUF 
GO  TO  445 

444  WRTTE(3,33)FBUF 

445  CONTINUE 
33  FORMAT(F) 

IF(INFOU).EQ.O)  GO  TO  60 
31  CONTINUE 

IF(IUNIT.EQ.3)GO  TO  888 
CLOSE(UNIT=2,DISP=’SAVE') 

GO  TO  889 

888  CLOSE(UNrr=3,DISP='SAVE’) 
c 

C  WHEN  ALL  DATA  ARE  ACQUIRED,  PRINT  THE  RAW  DATA  AND 

C  END. 

C 

889  CONTINUE 
TYPE  1026,INFO(1) 

1026  FORMAT!/’  A/D  SWEEP  ENDING  CODE  =  T3) 


CROSS.  FOR 

CROSSCORRELATION  SEQUENCE 
IN  ORDER  TO  CROSSCORRELATE  TWO  SIGNALS.TWO 
PROGRAMS  MUST  BE  EXECUTED  IN  THE  ORDER  THAT  THEY 
APPEAR  BELOW: 

(D...RUN  CROS1 
(2).. .RUN  CROSS 

THE  FIRST.’CROSl’ .CONVERTS  A/D.DATA  IS  GATHERED  FROM 
THE  FIRST  TWO  CHANNELS  ON  THE  A/D  BOARD  AND  8,000 
POINTS  ARE  STORED  IN  A  SINGLE  FILE  CALLED  FILE.DAT. 

THE  A/D  BOARD  SAMPLES  AT  A  RATE  OF  1MSEC, THEREFORE 
DATA  MUST  BE  GATHERED  FOR  8  SECONDS. 

THE  2ND, 'CROSS’. IS  THE  ACTUAL  CROSSCORRELATION. 

TWO  SIGNAL  FILES  ARE  CREATED  FROM  THE  ONE  CREATED 
BY  CROS1. 

THE  TWO  SIGNAL  FILES  ARE  MULTIPLIED  AND  THEN 
INTEGRATED  EACH  TIME  THE  PROGRAM  GOES  THROUGH  ITS 
MAIN  LOOP  ONE  VALUE  OF  THE  CROSS-CORRELATION 
COEFFICIENT  IS  STORED  IN.  THE  VARIABLE  LOCATION 
SUMT.EACH  POINT  CORRESPONDS  TO  A  DELAY  TIME  EQUAL 
TO  1MSEC.EXCEPT  FOR  THE  FIRST  WHICH  CORRESPONDS  TO 
A  DELAY  TIME  EQUAL  TO  0.0  MSEC  300  VALUES  OF  THE 
CROSS-CORRELATION  COEFFICIENT  ARE  PRINTED.THEREFORE 
A  TOTAL  DELAY  TIME  OF  299  MSEC  IS  SHOWN.THE  PROGRAM 
SEARCHES  THE  300  POINTS  AND  LOCATES  THE  MAXIMUM 
CROSS-CORRELATION  VALUE  AND  THE  TIME  CORRESPONDING 
TO  THIS  MAXIMUM. THE  PROGRAM  THEN  CALCULATES  A 
PARTICLE  VELOCITY  BASED  ON  A  2.16  FT  PROBE  SEPARATION. 
THIS  PROGRAM  ALSO  CALCULATES  THE  VOIDAGE(E),THE  GAS 
VELOCITY(VG).THE  SLIP  VELOCITY(VS),THE  ACTUAL  GAS 
RATE(SCF M A) ,THE  SOLID  FLOW  RATE(WS),AND  THE 
THE  LOADING. 

DURING  EXECUTION  A  NUMBER  OF  QUESTIONS  ARE  ASKED. 


DIMENSION  FBUF(2500),FBUF1(2500),RMULT(2200), 
1  SUMT(300) 

TYPE  2990 

FORMAT(5XTHIS  MONTH  IS?  (1-12)'$) 
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ACCEPT  2991.KMON 

2991  FORMAT(I) 

TYPE  2992 

2992  FORMAT(5X, TODAYS  DATE  IS?’S) 

ACCEPT  2993.KDATE 

2993  FORMAT(I) 

TYPE  2994 

2994  FORMAT(5X,THIS  YEAR  IS7’S) 

ACCEPT  299S.KYEAR 

2995  FORMAT(I) 

TYPE  29% 

29%  FORMAT(5X,  'PRINT  TRIAL  NUMBER’S) 

ACCEPT  3001.KTN 

3001  FORMAT(I) 

TYPE  3002 

3002  FORMAT(5X,’PRINT  TAPE  NUMBER’S) 

ACCEPT  3003.KTAPN 

3003  FORMAT(I) 

TYPE  3004 

3004  FORMAT(5X,TAPE  FOOTAGE.BEGINNING  NUMBER  ?’S) 
ACCEPT  3005.KBTAF 

3005  FORMAT(I) 

TYPE  3006 

3006  FORMAT(5X,’TAPE  FOOTAGE.END?’S> 

ACCEPT  3007,  KENT  A 

3007  FORMAT(I) 

TYPE  3008 

3008  FORMAT(5X.’TAPE  SIDE  7*S) 

ACCEPT  3009.KTS 

3009  FORMAT(I) 

TYPE  4000 

4000  FORMAT(5X, ’SOLID  FLOW  RATE  (LBS/ MIN)  ?’S) 
ACCEPT  4001, WS 

4001  *  FORMAT(F) 

TYPE  4002 

4002  FORMAT(5X,’AIR  HUMIDITY  ?’S) 

ACCEPT  4003, AH 

4003  FORMAT(F) 

TYPE  3010 

3010  FORMAT(5X,’WHAT  IS  THE  TURBINEMETER  READING 
1  (SCFM)  7’S) 

ACCEPT  301LSCFM 

3011  FORMAT(F) 

DWELL=.001 
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TYPE  2003 

2003  FORMAT(5X,’SOLID  DENSITY  (LBS/FT-3)  ?*$) 

ACCEPT  2004.DEN 

2004  FORMAT(F) 

. . . . . . * . . 

C  THIS  STATEMENT  IS  USED  TO  CALCULATE  THE  ACTUAL  GAS 
C  FLOW  RATE.SCFM  IS  THE  TURBINEMETER  READING. 

C 

SCFMA=SCFM 

C 

. . . . •••*•••• . . . . . 

c 

. . . . . . . . 

C  THIS  STATEMENT  IS  USED  TO  CALCULATE  THE  ACTUAL 
C  LOADING  (ALOAD). 

C 

ALOAD=379*WS/(SCFMA*29) 

C 

c 

c 

. . * . . . 

c 

. . . . . . . . 

C  THE  OPEN  STATEMENT  READ  IN  THE  TWO  SIGNAL  FILES. FBUF 

C  IS  THE  VARIABLE  NAME  FOR  SIG1.DAT.AND  FBUF1  IS  THE 

C  VARIABLE  NAME  FOR  SIG2.DAT. 

C 

OPEN(UNIT=2.NAME=,FILE.DAT\TYPE=,OLD\DISP=,SAVE’) 

DO  6  ID= 1,2500 
READ(2.7)FBUF(ID) 

READ(2,7)FBUF  KID) 

6  CONTINUE 

41  PRINT  42 

42  FORMATCRING  PROBES  WERE  USED’,/) 

46  ADD1=0.0 

ADD2=0.0 

DO  45  1=1,2500 

ADD1=ADD1+FBUF(I) 

ADD2=ADD2+FBUF1(I) 

45  CONTINUE 

IF(ADDl.GE.0.0)GO  TO  47 
DO  48  1=1,2500 
FBUF(I)=-FBUF(I) 

48  CONTINUE 
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47 


145 


49 


43 

O 

C 

c« 


90 


100 

80 


c 

c« 

c 

c* 

c 

c 


IF(ADD2.GE.O.O)GO  TO  49 
DO  145  1=1.2500 
FBUF1(I)=-FBUF1(I) 

CONTINUE 

ADD1=0.0 

ADD2=0.0 

DO  43  1=1.2500 

IF(FBUF(I).LE.0.0)FBUF(I)=0.0 

IF(FBUF1(I).LE.0.0)FBUF1(I)=0.0 

CONTINUE 


DO  80  1=1,625 
TOTAL1=0.0 
TOTAL2=0.0 
DO  90  J=l,2 
M=J+(I-1)*2 

TOTALl=TOTALl+FBUF(M) 

TOTAL2=TOTAL2+FBUFl(M) 

CONTINUE 

AVG1=TOTAL1/2.0 

AVG2=TOTAL2/2.0 

DO  100  K=l,2 

M1=K+(I~1)*2 

FBUF(M1)=AVG1 

FBUF1(M1)=AVG2 

ADD1=ADD1+FBUF(M1) 

ADD2=ADD2+FBUF1(M1) 

CONTINUE 

CONTINUE 

THRESl=0.5*ADDl/2500.0 
THRES2=0. 5  »ADD2/ 2500.0 


THE  FOLLOWING  STATEMENT  SET  THE  THRESHOLD  VALUE. 
DO  10  1=1.2500 

IF(FBUF1(I).LE.THRES2)FBUF1(I)=0.0 

IF(FBUF(I).LE.THRES1)FBUF(I)=0.0 

CONTINUE 


10 

C 


. . * . 

C  THIS  IS  WHERE  THE  PROGRAM  MULTIPLIES  AND  INTEGRATES 
C  THE  TWO  SIGNAL  FILES. 

C 

DO  11  J=1.300 
DO  12  1=1.2200 
K=I+J 

RMULT(I)=FBUF(K)*FBUF1(I) 

12  CONTINUE 

SUM =0.0 
DO  60  1=1.2199 

C12=(RMULT(MhRMULT(I))*DWELL/2.0 
SUM=SUM+C12 
60  CONTINUE 

SUMT(J)=SUM 
11  CONTINUE 

C 

. . . . . . . . . . . 

c 

. . . . . . 

C  THIS  STATEMENT  TAKES  THE  GAS  FLOW  RATE  AND  DIVIDES 

C  IT  BY  THE  PRODUCT  OF  THE  CROSS  SECTIONAL  AREA  AND 

C  THE  REQUIRED  TIME  CONVERSION  FROM  MINUTES  TO 

C  SECONDS. 

C 

GV=SCFMA/1.3712427 

C 

. . . . . . . . . . . . . 

c 

. . MM... . . . . . . 

PRINT  101.KMON.KDATE.KYEAR 

101  FORMAT(/.’DATE’,2X,I,2X.I.2X,I,/  /) 

PRINT  102.KTN 

102  FORMATCTRIAL  NUMBER’, 2X.I.//) 

PRINT  103.KTAPN 

103  FORMATCTAPE  NUMBER’, 2X.I,//) 

PRINT  104.KBTAF.KENTA 

104  FORMATCTAPE  FOOTAGE'.lX,I.’-\I.//) 

PRINT  105.KTS 

105  FORMATCTAPE  SIDE\2X,I.//> 

PRINT  106 

106  FORMAT(50X,’CROSSCORRELATION  DATA'.//) 


369 

C 

o» 

c 

c» 

c 

c 

c 


831 


832 

833 

134 


500 

C 

o* 

c 

o* 

c 

c 

c 

900 


350 

368 

C 


PRINT  369.SUMT 
FORMAT(10(F10.7,2X)./) 


THE  FOLLOWING  STATEMENTS  SEARCH  THE  CROSS¬ 
CORRELATION  DATA  TO  LOCATE  TIMES  WHEN  MAXIMUMS 
OCCUR. 

TAU=2166.67/GV 
KOUNT=0 
TIME=0.0 
DO  500  1=1.298 
TIME=TIMEfl.0 

IF(SUMT(I).GRSUMT(I+1))  GO  TO  134 
IF(SUMT(I).LT.SUMT(I+1))  GO  TO  831 
CONTINUE 

IF(SUMT(M).LT.SUMT(I+2))  GO  TO  134 
IF(SUMT{I+l).GT.SUMT(I+2))  GO  TO  832 
PRINT  833,TIME,SUMT(M) 

FORMATCA  MAXIMUM  OCCURED,TIME(MSEC.)=’.F5.1.5X, 
F11.7./) 

IF(TIME-LT.TAU)  GO  TO  500 
IF(KOUNT.EQ.l)  GO  TO  500 
N=I 

KOUNT=KOUNT+l 

T=TIME 

CONTINUE 


THESE  STATEMENTS  LOCATE  THE  MAXIMUM  CROSSCORRELATION 
VALUE  AND  THE  CORRESPONDING  DELAY  INTERVAL. 

PMAX=SUMT(N) 

KDT=1 

DO  368  I=N,300 

IF(SUMT(I).LT.PMAX)  GO  TO  350 
PMAX=SUMT(I) 

TIMAX=T 

T=T+KDT 

CONTINUE 


C***** . *••• . * . . . . ****** . 

C  THIS  STATEMENT  CALCULATES  THE  PARTICLE  VELOCITY. 

C  SINCE  TIME  IS  IN  SECONDS  INSTEAD  OF  MSEC.  A  FACTOR 
C  OF  1,000  APPEARS  IN  THE  NUMERATOR. 

C 

PV =2167.67  /TIM  AX 
C 

. . . . ••** 

c 

. . . . . . . . * 

C  THE  NEXT  TWO  STATEMENTS  ARE  USED  TO  CALCULATE  THE 
C  VOIDAGE(E). 

C 

EA=0.72963*ALOAD*SCFMA/(PV*DEN) 

E=1-EA 

C 

. . . . •••••*••••••* 

c 

. . ****•••• . MM. . 

C  THIS  STATEMENT  CALCULATES  THE  SLIP  VELOCITY. 

C 

SV=(GV/E)-PV 

c 

. . * . . . . . . 

c 

. . * . * . * . 

PRINT  132.  AH 

132  FORMATOAIR  HUMIDITY(%)\2X,F5.2.//) 

PRINT  131.ALOAD 

131  FORMAT(’LOADING’,2X,F6.3,/  /) 

PRINT  358.SCFMA 

358  FORMAT!’ AIR  FLOW  RATE(SCFM)’,2X,F6.3,//) 

PRINT  701.E 

701  FORMAT(’VOIDAGE’,2X,’E=’,F,/  /) 

PRINT  702, WS 

702  FORMATCSOLID  FLOW  RATE  (LBS/MIN)\2X.’WS=’,F5.3.//) 

PRINT  703.GV 

703  FORMATCGAS  VELOCITY  (FT/SEC)\2X,*VG=\F6.3,//) 

PRINT  704, PV 

704  FORMATCPARTICLE  VELOCITY  (FT/SEC)’,2X,’VP=\ 

1  F6.3,//) 

PRINT  705.SV 
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705  FORMATCSLIP  VELOCITY  (FT/SEC),,2X,,VS=,,F6.3.//) 
CLOSE(UNIT =3) 

CLOSE(UNIT=4) 

STOP 

END 
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C  PROGRAM  TO  CALCULATE  CHOKING  AND  SALTATION 
C  VELOCITIES 
C 

100  PRINT*, ’INPUT  THE  PARTICLE  DENSITY* 

READvRP 

RF=1.2 

XMU=l.E-5 

G=9.8 

D=0508 

A=3.1415*D**2/4 

PRINT*. 'INPUT  THE  PARTICLE  DIAMETER’ 

READ*. DP 
DP=DP*l.E-6 

PRINT*. 'INPUT  THE  MASS  FLOW  RATE’ 

READ*,WS 

FPC=6.81E5*(RF/RP)**2.2 

PRINT ‘.’INPUT  AN  INITIAL  GUESS  FOR  UFC 

READ*, GUESS 

UFL=GUESS 

DO  10  1=1.500 

UT=.153*G**.71*DP**1.14*(RP-RF)**.7l 

UT=UT/(RF**.29*XMU**.43) 

UF=(2*G*D*((l-(WS/((UFL-UT)*RP*A)))**-4.7-l)/FPC)**.5 

UF=UF+UT 

ERR=ABS(UF-UFL) 

IF(ERR.LT..0001)  THEN 

GOTO  20 

ENDIF 

UFL=UF 

10  CONTINUE 

PRINT*, 'DID  NOT  CONVERGE  FOR  YANG’ 

GOTO  30 

20  PRINT*, ’UFC  YANG  =  ’,UF 

30  CONTINUE 

DP=DP*1000. 

UGSALT=WS/(RF*A)*10.**(1.44*DP*-1.96) 

UGSALT=UGSALT*(G*D)**((1.18*DP+2.5)/2) 

UGSALT=UGSALT**(l/(l.l*DP+3.5)) 

PRINT*, ’SALTATION  VELOCITY  BY  RIZK=  ’.UGSALT 
DP=DP/1000. 

PRINT*, ’INPUT  AN  INITIAL  GUESS  FOR  OWENS  SALTATION’ 

READ*,UFLO 

DO  40  1=1,500 

FG=.0014+.125/((RF*D*UFLO)/XMU)**,32 


UGSALT=(.01*(2*RP*G*DP)/(RF*FG))**.5 

ERR=ABS(UGSALT-UFLO) 

IF(ERR.LT.,0001)  THEN 

GOTO  60 

ENDIF 

UFLO= U  GS  ALT 
CONTINUE 

PRINT*, ’DID  NOT  CONVERGE  FOR  OWENS’ 

GOTO  50 

PRINT*, ’SALTATION  VELOCITY  BY  OWENS  =  ’.UGSALT 
CONTINUE 

PRINT*,’  INPUT  AN  INITIAL  GUESS  FOR  ROSE’ 

R£AD*,UFL 
DO  70  1=1.500 
WG=RF*A*UFL 

UGCR=3.2*(WS/WG)**.2*(D/DP)**.6*(RP/RF)*— .7 

UGCR=UGCR*(UFL/(G*D)**.5)**.5 

UGCR=UGCR*UT 

ERR=ABS(UGCR-UFL) 

IF(ERR.LT..0001)  THEN 

GOTO  80 

ENDIF 

UFL=UGCR 

CONTINUE 

PRINT*,’DID  NOT  CONVERGE  FOR  ROSE  AND  DUCKWORTH’ 

GOTO  90 

PRINT*. ’CHOKING  VLOCITY  BY  ROSE  AND  DUCKWORTH  =  ’.UGCR 
CONTINUE 

PRINT*, 'INPUT  A  1  TO  CONTINUE’ 

READ*,Z 

IF(Z-EQ.l.)  THEN 

GOTO  100 

ENDIF 

STOP 

END 


LINEAR  STABILITY 

DT=PIPE  DIAMETER 

DP=PARTICLE  DIAMETER 

UPS=PARTICLE  VELOCITY 

UFS=GAS  VELOCITY 

RP=PARTICLE  DENSITY 

RF=FLUID  DENSITY 

PRESS = PRESSURE  DROP 

XM1.XM2  ARE  THE  EIGENVALUES 

IRUN  IS  THE  EXPERIMENTAL  RUN  NUMBER 

OPEN(UNIT=10.FILE=’INPUT.DAT’) 

DT=0.0508 

RF=1.2 

GRAV=9.8 

PRINT*, 'INPUT  PARTICLE  DIAMETER  IN  MICRONS’ 

READ*, DP 
DP=DP*l.E-6 

PRINT*, 'INPUT  ANGLE  OF  INCLINATION  IN  DEGREES’ 

RE  AD*.  ANGLE 

THETA= ANGLE/ 180.  *3. 1415 

PRINT*,SIN(THETA),'  IS  THE  SIN  OF  THETA’ 

GET  CONSTANTS  FOR  PARTICLE  VELOCITY  AS  FUNCTION 
OF  GAS  VELOCITY  (LINEAR) 

PRINT*, ’INPUT  THE  SLOPE’ 

READ*  .SLOPE 

PRINT*, ’INPUT  THE  INTERCEPT’ 

READ*.  INTER 

OPEN(UNrT=ll,FILE=’OUTPUT.DAT) 

WRITE(11,200) 

WRTTE(11,210) 

WRITE(11,220) 

WRITE(11,210) 

WRITE(11,200) 

WRJTE(ll,2lO) 

FORMAT(5(’ _ ’)) 

FORMATC  ’) 

FORMAT(4X,’RUN’,7X,’Ug’,8X,’Uf’,10X.’Ml’,10X,’M2’) 

READ(10,110)NDAT,RP,IRUN 

FORMAT(I,F,I) 

IF(NDAT.EQ.999)  THEN 

GOTO  1000 

ENDIF 

DO  10  I=1,NDAT 
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READ<10.100)WS.UFS,UPS.PRE,DUMMY1,DUMMY2.DUMMY3 
100  FORMAK7F) 

RE=DT*UFS*3.28*0.07476/(1.24E-5) 

FG=0.0791/RE**0.25 
DPG=2*FG*UFS**2*0.07476»10/0. 170583 
DPG=DPG*760  /  (32. 174*  14. 696*  144) 

D*FF=DUMMY2-DPG 

PRESS=PR£-DIFF 

C  CONVERSION  TO  SI  UNITS 

WS=WS*7.56E-3 
UFS=UFS*0.3048 
UPS=UPS*0.3048 
PRESS=PRESS*43. 74097 
RP=RP*  16.02 
UPS=SLOPE*UFS+ INTER 
E-1-(4*WS/(3.1415*DT**2*UPS*RP)) 

C  CALC  DRAG  COEFF 

REP=RF»DP*(UFS/E-UPS)/l.E-5 

CDS=18.5*E**-4.7/REP**.6 

UFS=UFS/E 

FS=0.0285*(GRAV*DT)**.5/UPS 
C  GET  CONSTANTS  FOR  LINEAR  ANALYSIS 

A=3*CDS*RF/(4*RP*DP) 

B=GRAV  *SIN(THETA)+PRESS/RP 
C~2*FS/DT 

D '  JrRAV*SIN(THETA)+PRESS/RF 
EE=2*FG/DT 

A0=A*(UFS-UPS)**2-B-C*UPS**2 

A1=2*A*(UFS-UPS) 

A2=-2*C*UPS 

B0=-A*(UFS-UPS)**2-D-EE*UFS*»2 

B1=-2»A*(UFS-UPS) 

B2=-2*EE*UFS 

XL1=(A1-B1-B2-A2) 

XL2=A2*B1+A2*B2-A1*B2 

XL3=A1*B0-A0*B1-A0*B2+A2*B1*UPS-A1*B2*UPS+A2*B2*UPS 
C  D2+L1XD1+L2=L3 

XM1=-XL1+(XL1**2-4*XL2)**.5 
XMl=XMl/2 

XM2=-XL1-(XL1**2-4*XL2)**.5 
XM2=XM2/2 
UFS=UFS*E 

WRJTE(11,Z50)IRUN,UFS,UPS,XM1.XM2 
FORMAT(4XJ3,5X,F5.2,5X.F5.2,5X,F8.3,5X,F8.3) 
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IRUN=IRUN*1 
RP=RP/ 16.02 
10  CONTINUE 

WRITE(1U60) 
260  FORMAT(’l’) 

GOTO  5 
1000  STOP 


PROGRAM  TO  CALCULATE  PRESSURE  DROP 
RATIOS  WITH  YANG  AND  KONNO  SAITO 
DT=PIPE  DIAMETER 
DP=PARTICLE  DIAMETER 
UPS=PARTICLE  VELOCITY 
UFS=GAS  VELOCITY 
RP=  PARTICLE  DENSITY 
RF=FLUID  DENSITY 
PRESS=PRESSURE  DROP 
IRUN  IS  THE  EXPERIMENTAL  RUN  NUMBER 
DPKS  =  PRESSURE  DROP  BY  KONNO-SAITO 
DPY  =  PRESSURE  DROP  BY  YANG 

OPEN  INPUT  FILE 

OPEN(UNIT=10.FILE=TNPUT.DAT) 

DT=0.0508 

RF=1.2 

GRAV=9.8 

GET  DP  AND  ANGLE  FROM  SCREEN 

PRINT*. 'INPUT  PARTICLE  DIAMETER  IN  MICRONS’ 

READ*. DP 

DP=DP*l.E-6 

PRINT • , ’ INPUT  ANGLE  OF  INCLINATION  IN  DEGREES’ 

READ*  .ANGLE 

THETA=ANGLE/180.*3.1415 

OPEN  OUTPUT  FILES 

OPEN(UNIT=ll,FILE=’ROUT.DAT) 

OPEN(UNIT=12.FILE=’ROUTY.DAT’) 

OPEN  (UNIT = 13,F  ILE= ’  MEAN  .DAT’  ,ACCESS=’  APPEND’ ) 
FORMAT!’  ’) 

READ  DATA  AND  CHEK  FOR  END  OF  FILE 

R£AD(10,110)NDAT,RP.IRUN 

PRINT*. ’ND AT  IS  ’.NDAT 

FORMAT(I.F.I) 

IF(NDAT.EQ.999)  THEN 

GOTO  1000 

ENDIF 

DO  10  I=1,NDAT 

READ(10,100)WS,UFS,UPS,PRE.DUMMY1,DUMMY2,DUMMY3 

FORMAT(7F) 

RE=DT*UFS*3.28*0.07476/(1.24E-5) 

FG=0.0791/RE**0.25 
DPG=2*FG*UFS**2*0.07476*  10/0. 170583 


'■m  ^ 
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DPG=DPG*760/(32. 174*14.696*144) 

DIFF=DUMMY2-DPG 
PRESS=PRE-DIFF 
CONVERSION  TO  SI  UNITS 
WS=WS*7.56E-3 
UFS=UFS*0. 3046 
UPS=UPS*0.3048 
PRESS=PR£SS*43.73097 
RP=RP*  16.02 

E=1-(4*WS/(3.1415*DT**2*UPS*RP)) 

UFS=UFS/E 

CALC  TERMINAL  VELOCITY 
UT=0.153*GRAV**.71*DP**1.14*(RP-RF)**.71 
CALC  PRESSURE  DROP  BY  KONNO  SAITO  AND 
GET  THE  RATIO  RATIOK 
UPKS=E*UFS-UT 

FSKS=0.0285*(GRA  V  •  DT)*  *.5  /  UPKS 

EKS=1-(4*WS/(3.1415*DT**2*(RP-RF)*UPKS)) 

DPKS=(1-EKS)*GRAV*SIN(THETA)*RP+RF*EKS*GRAV*SIN(THETA) 

DPKS=DPKS+2*FG*RF  *UFS**2/DT 

DPKS=DPKS+2*FSKS*RP*(1-EKS)*UPKS**2/DT 

RATIOK=PRESS  /  DPKS 

DEVKS=ABS((PRESS-DPKS)  /PRESS)*  100 

DECIDE  WHICH  YANG  CORRELATION  TO  USE 

VERTICAL  OR  HORIZONTAL 

IF(ANGLE.LT.50.)  THEN 

GOTO  200 

ENDIF 

CALC  PRESSURE  DROP  BY  YANG 
FOR  VERTICAL  PIPE  AND 
PRINT  RESULTS  TO  THE  OUTPUT 
FILE 
UPY=1. 

EY=1-4*WSV(3.1415*DT**2*(RP-RF)*UPY) 

FSY=((l-EY)*UT/(UFS-UPY))**-0.979 

FSY=FSY*(1-EY)*.00315/EY**3 

UPYC=UFS-UT*SQRT(((1+FSY*UPY**2/(2*GRAV*DT))*EY**4.7)) 

ERROR=ABS(UPYC-UPY) 

IF(ERRORLT.O.OOOl)  THEN 

GOTO  180 

ENDIF 

UPY=UPYC 

GOTO  150 

DPY =(1-EY)*GRA  V*RP+RF*EKS*GRA  V 
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DPY=DPY+2«FG*RF*UF**2/DT 
DPY=DPY+2»FSY»RP*(1-EY)»UPY**2/DT 
RATIO  Y =PRESS  /  DP  Y 
DE  VY= ABSUPRESS-DP  Y)  /  PRESS)  •  100 
RP=RP/ 16.02 

WRITE(11.350)UFS,RATIOK. 

WRJTE(12.350)UFS.RATIOY 
WRITE!  13.370)DE  VKS.DE  VY 
GOTO  10 

C  CALC  PRESSURE  DROP  BY  YANG  FOR 

C  THE  HORIZONTAL  PIPE  AND  PRINT 

C  RESULTS  TO  THE  OUTPUT  FILES 

200  UFY=1. 

205  EY=1-4»WS/(3.1415*DT*»2*(RP-RF)*UPY) 

IF(EY.GT.0.9999)  THEN 

EY=.  99999 

ENDIF 

FSY=.0293*(1-EY)/EY»*3 
FSY=FSY*((1-EY)»UFS/(GRAV«DT)**.5)*— 1.15 
UPYC=UFS-UT«SQRT(((1->FSY»UPY»*2/(2*GRAV*DT))»EY**4.7)) 
IF(UPYGLT.O.O)  THEN 
UPYC=ABS(UPYC) 

ENDIF 

ERROR=ABS(UPYC-UPY) 

IF(ERROR-LT.O.OOOl)  THEN 

GOTO  300 

ENDIF 

UPY=UPYC 

GOTO  205 

300  DPY=2*FG*RF*UFS*»2/DT 

DPY=DPY+2»FSY«RP»(1-EY)»UPY»»2/DT 

RATIO  Y =PRESS  /  DP  Y 

DEVY= ABS!  (PRESS-DPY)  /  PRESS)  •  100 

UFS=UFS*E 

WRITE(11.350)UFS.RATIOK 
WRITE(12,350)UFS.RATIOY 
WRITE!  13.370)DE  VKS.DE  VY 
350  FORMAT(F.’,’,F) 

15  RP=RP/ 16.02 

10  CONTINUE 

WRJTE(11,360)IRUN 
WRITE!  12.360)IRUN 
WRITE!  13. 360)IRUN 

360  FORMATCRUNS  ABOVE  BEGAN  AT  RUN  M) 
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